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I. INTRODUCTION 


This research has grown out of previous work (Craig, 1982a) In Which 
It was shown that landforms have a characteristic scale which can be 
determined by analysis of the degree of relation between adjacent slopes. 
Later (Craig 1982b) it was shown that the existence of this relation 
between adjacent slopes could be understood in light of certain "classical'' 
theories of slope process and form which had been expressed as differential 
equations. Moreover, it now appears possible to infe*' process rates 
from observed forms in a fairly objective and precise manner. 

Such work relates to remote sensing problems of interest to NASA 
because it has previously been shown that the above mentioned 'scale' of 
the landform impacts reflectance from that landform so as to yield a 
characteristic scale and dependence between adjacent pixels (picture 
elements) of LANDSAT data (Craig, 1979, 1981; Craig and Labovitz, 1980). 

An even more intimate relation exists between the geobotanical 
studies of Dr. Mark Labovitz of NASA's Goddard Space Flight Center and 
the geomorphic theories outlined above. In these studies it is desired 
to find ways to recognize the surface expression of significant ore bodies 
using remotely sensed data. Emphasis is now being placed on the eastern 
United States. Here the ore body is likely to be obscured by regolith 
and vegetation. Thus one important aspect of the study must be to under- 
stand to what extent processes at work on the regolith are likely to 
influence the rate and mode of movement of materials from the ore body. 

This study addresses five principal objectives related to the fore- 
going points. The first objective is to study and to characterize the 
topographic complexity of each of eight physiographic provinces in the 
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•astern half of the United States. This study is United to eight provinces 
because of its exploratory nature. These particular provinces (see Figure 1) 
were chosen because they are of greatest Interest at present In the ongoing 
geobotanical invest igatievs. They are listed In Table I and the boundaries 
follow those of Fanneman (1938) . Topographic complexity Is here defined 
as the degree of relation between contiguous elope segments In a traverse 
as measured by the autocorrelation vBox and Jenkins, 1970). An important 
aspect of this of the work is to determine an appropriate sampling 

schem,-.' for each *<',te. Thi;- is discussed in the section on methods (page 33). 

Once '.o* toy„grf*?u.»c complexity has been characterized at each site 
the refund step is ~o corny re the variability of autocorrelation within a 
oroall area (here defined as a single 74’ quadrangle) to the variability 
at widely separated and diverse areas within the same physiographic 
region. The objective of this step is to obtain some measure of the 
degree of uniformity of the autocorrelation (and hence presumably of the 
processes) which can be expected to be encountered within a given physio- 
graphic province. Since each province is supposed to represent a homo- 
geneous structural, lithologic, climatic and geomorphic system it is 
reasonable to expect that the measures obtained for topographic complexity 
will be uniform within each province. If this can be shown to he so it 
will simplify the design of the geobotanical studies. 

Whereas uniformity within provinces might be expected it is quite 
likely that provinces will differ in the extent to which each process is 
active both in an absolute and a relative sense (for a discussion of the 
processes being detected see Chapter II, "Fundamental Theory"). For 
example, we can expect slope wash to be more important in a province 
with a climate typified by Intense rain showers and with a thin vegetative 
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Tabic I. Physiographic province* actually studied. 

1. Interior Low Plateaus 

2. New England province 

3. Piedmont province 

4. Blue Ridge province 

5 . Ozark Plateaus 

6. Valley and Ridge 

7 . Appalachian Plateaus 

8. Ouachita province 
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cover. On the other hind, creep night be aore prevalent where the regolith 
ie thick and teapereture extremes aore severe. Thus, the third aajor 
objective of this study is to compare and contrast the variability of 
autocorrelation across the eight physiographic regions. The results of 
this investigation should provide a guide to the extent of effort required 
in the geobotanlcal studies in each province. 

Even if it is demonstrated that each physiographic region is homo* 
geneous in topographic complexity within its boundaries and that there are 
significant differences between regions that does not imply that each 
province is distinct from every other one. It is not unlikely that we 
could find that provinc*- :: and Ye*, similar to one another but distinct 
from provinces Z and U (which also define a group). Thus, a fourth aajor 
goal is to partition the total study area into subareas homogeneous in 
terrain complexity (autocorrelation properties). This study will also 
allow for the possibility that region A of province X may be more like 
region C of province Y than it is like region B of province X. Once 
these homogeneous subareas have been defined (if any exist at all) it 
will be easier to specify which geobotanlcal sampling procedure is 
appropriate in a given area. 

The last point Implies that there are distinct strategies appropriate 
for areas depending upon the “terrain complexity" of that area. Thus, an 
appropriate question is to what extent should the differences in terrain 
complexity be translated into different geobotanlcal sampling techniques. 
Therefore, the fifth and final goal of this research is to show the rela- 
tion between the complexity measured, the geomorphic process mix implied 
and the way in which geobotanlcal information will be modified into a new 
and more or less recognizable entity. 


* 
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A discussion of the underlying theoretical framework of this study 
in the next section will be followed by a statement of the sampling method 
employed. Results are given next and are shown to be remsrkably consistent 
with general expectations and particularly helpful in understanding 
geomorphic processes. This is followed by "Interpretation" both geomorphic 
and geobotanlcal. A section summarizing the major conclusions of this 
study (Chapter VI), is followed by some recommendations for future work 
(Chapter VII). 


9 


II. FUNDAMENTAL THFO RY 

The elope forming processes of interest are those which are dominantly 
transport limited (Carson and Kirkby, 1972) since these are most expected 
in regions typified by significant colluvial mantles as is common in much 
of the eastern U.S. Such processes are here divided into two groups and 
called, for brevity, slope wash and creep. These terms are shorthand 
for process groups which include solution on the one hand and solifluction 
on the other. 

Such groups are distinguished by the underlying controls upon their 
rate of activity. Both are driven by gravity but the surface wash group 
is assumed to act at a rate proportional to the angle of slope. Steeper 
slopes will tend to display greater activity. On the other hand rates of 
creep will be influenced by the rate of curvature of the slope. It will 
be most active where the difference in slope angles between adjacent slope 
facets is the greatest. 

These ideas can be made more intuitive when the characteristics of 
the processes are considered. Slope wash processes take material from 
the slope in a direct way. Once the material has been included in the 
transporting agent (usually water) it usually remains there until that 
agent has debouched from the slope. Thus, it is sufficient to set the 
transporting agent in motion. On the other hand creep acts in an incremental 
and cumulative fashion. Removal of material from one slope facet does 
not embed it in an agent that will totally remove it from the slope. 

Rather that material moves directly to the adjacent slope facet, which is 
now burdened with this additional load to be removed plus any regolith to 
be removed at that spot. Thus, the cumulative nature of creep processes 
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is such that creep can only be active on the vhole slope when the slope 
angles are continuously increasing. The slope profile must be convex 
outward . 

These theoretical notions have been stated more succinctly as a 
set of differential equations (Culling, 1965). These equations relate the 
rate of change of elevation at each point of the slope to the relevant 
slope parameter through a constant of proportionality. For slope wash 
the form is: 

iz „ b . h. 

3t 3x (1) 

and for creep it is: 



In both of these equations x is taken to be the distance from the drainage 
divide of the slope, y is elevation and t is time. As can be seen these 
equations represent the development of a two-dimensional slope profile. 

They also are exact in that they fail to consider the possible influence 

of some other agent upon the slope profile. Such an influence could be 

most easily allowed by the addition of a random term to the original 
equation. 

Another drawback to the existing equations is the lack of an objec- 
tive means of estimating the values of the coefficients a and b. Typically 
if they are desired the equations must be solved for a variety of values 

and the correct value chosen by visual comparison to the slope of Interest. 

Although not stated explicitly it appears to be reasonable to expect both 
values to be positive and less than one. 

The separate equations have been combined into a single equation 
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(Hirano, 1975) which, under what appear to be reasonable assumptions, can 
be simplified to: 

3y 3y 3 2 y 

— “ b — + a — x (3) 

3t 3x 3x 

The terms carry the same meaning as previously. This equation predicts 
the combined effect upon the slope profile of the two distinct processes. 
Again there is no formal means to include the effects of other agents, 
and there is no objective method of estimating a or b. However, these 
equations, most importantly the last one, can be converted to a discrete 
form and by making use of the ergodic hypothesis (substitution of space 
for time) can be related to existing profiles at specific locations. 
Moreover, these discrete forms have the advantage that they are members 
of a family of discrete equations which has been intensely studied in 
recent years (Box and Jenkins, 1970) and for which extremely sophisticated 
and comprehensive statistical tests are available. 

The discrete form is achieved in this manner. Take 


3y 

3y 

a 2 y 


— 

■ b — + 

• — ; 

(4) 

St 

3x 

3x 2 



which is equivalent to: 

E*(t) * bE*(x-l) + a E 2 (x-1) (5) 

and, under the ergodic hypothesis 

E 1 (t) - E*(x) 


( 6 ) 


* &4I 
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so that 

E 1 (x) - b E 1 (x-1) + a E 2 (x-1) 


(7) 


which is equivalent to 

a{(E(x-l) - E(x-2)J - [E(x-2) - E(x-$J} + bfE(x-l) - E(x-2)] 

- (b+a)[E(x-l) - E(x-2)] - a[E(x-2) - E(x-3)J 

» (b+a) E* (x-1) - a E 1 (x-2) (8) 

and if we let 


■ a + b 
$2 * “• 


(9) 


and add a tern A(X) representing random effects introduced at point X 
independent of the agents being modelled we obtain an equation which is 
a member of the family of AutoRegressive-Integrated-Moving Average or 
ARIMA models. 


E d (x) ■ <J>jE d (x-l) + (x— 2 ) + . . . + 4ipE^(x-p) + 

A(X) - OjAtx-l) - © 2 A(x-2) - ... - 0 q A(x-q) 


( 10 ) 


where 


E d (x) - E d " 1 (x-1) - E d_1 (x-2) 


(ID 


The values of the parameters p, d and q define the order of the model. 
For the slopes the model is an ARIMA (2,1,0). 

There are a number of considerations of importance in the study of 
landforms using the Box and Jenkins (1970) methodology. Generally, the 
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study Is divided into three steps. The first is the identification of 
the appropriate model, that is the correct values of p, d and q for the 
series. It is suggested here that 

p - 2 

d « 1 (12) 

q ■ 0 

will be appropriate for almost any traverse of elevations observed in the 
eastern U.S. However, this assumption can be checked by using available 
identification procedures. 

Once the order of the model has been chosen the next step is to 
estimate the values of the component coefficients, in this case <J>j and 
Following this step comes diagnostic checking . In this stage a number of 
specific tests are available to determine if the model chosen adequately 
fits the empirical data. Each of these steps has been followed for the 
traverses measured in this study. Before these are reported a number of 
elements of these steps should be presented. 

Identification of the relevant model (i.e. the order of p, d and q) 
is most adequately done through the use of two items called the autocor- 
relation function and the partial autocorrelation function. Autocorrela- 
tion is the degree of relation between a sequence and itself at 
neighboring points at some fixed distance, k. Computing this correlation 
at successive values of k yields, a set of autocorrelations which can 
be plotted versus k to yield the autocorrelation function (ACF) . For 
an ARIMA (2,1,0) the ACF can take on any of four fundamental forms which 
vary according to the values of and These forms are shown in 


figure 2. 
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Figure 2. Examples of the four major different forma possible in an 

Autocorrelation Function of an ARIMA (2,1,0) model. Letters 
are keyed to areas of figure 5. 
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Other models with some other velue of p (i.e. p 4 2) would show 
distinctive forms of the ACF. The ACF of a model with p ■ 1 could be 
confused for those shown since it is equivalent to an AR1MA (2,1,0) with 
4*2 “ 0. However, each of these alternative possibilities can be detected 
through the use of the second item of diagnosis, the partial autocorrelation 
function. 

Partial autocorrelation is closely analogous to classical partial 
correlation. A simple example will show the value of its use. Suppose 
it is found that adjacent observations tend to be highly autocorrelated. 
Thus, point x influences point x + 1 while point x + 1 influences point 
x + 2. It is easy to see that point x will show a relation to point x + 2 
simply due to the carry-over effect through point x + 1 even if no actual 
two step autocorrelation exists. In order to test for true two step 
autocorrelation, it is first necessary to remove the effect of the one 
step autocorrelation carry-over. Such correction produces the partial 
autocorrelation at lag two. Similar corrections are done for each lag in 
order to remove the effects of autocorrelation at lower lags. The resulting 
values, when plotted versus the appropriate lag is called the partial 
autocorrelation function (PACF). 

Representative PACF's of ARIMA (2,1,0) models are shown in figure 3. 

In general, the PACF of an ARIMA (p,d,q) will have p-many significant 
values followed by values not significantly different from zero. It is 
the combination of these forms of the PACF together with an ACF as in 
figure 2 which is diagnostic of an ARIMA (2,d,0). Comparable diagnostic 
tools are available for other models and further details can be found in 
Box and Jenkins (1970). 

Following the diagnostic stage, estimates of the values of the $'s 
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Figure 3. Examples of the four major different forms possible in a 

Partial Autocorrelation Function of an ARIMA (2*1,0) model 
Letters are keyed to areas of figure 5. 
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and 0's arc obtained uaing a least squares method which yields an approxi- 
mate maximum likelihood solution. The technique is iterative and is most 
efficient if reasonable preliminary estimates are used as starting values. 
For the ARIMA (2,1,0) these can be obtained from the first two values of 
autocorrelation (that at lags 1 and 2, labelled Pj and P 2 using the 
following equations 


2 

♦j - P 1 (l-P 2 )/(1-Pi> 

2 2 
0 2 “ (P 2 "Pl^ ' (l“Pl) 


(13) 


Certain values of and $ 2 are not feasible in that they can only 
arise from a non-stationary process, that is one for which a mean value 
does not exist. Indeed, it appears that traverses of elevation do 
represent a non-stationary process; however, slopes, the series of 
interest in this study, are almost surely stationary. Three inequalities 
define the admissible region for the parameters and $2 which will 
yield stationary slope series 


4 2 + < 1 

* 2 - 4>i < 1 
- 1 < *2 < 1 

this means the parameters and $2 roust fall within the triangular region 
shown in figure 4. 

The four types of ACF and PACF shown in figures 2 and 3 arise when 
and $2 take on specific values. These can best be seen in figure 5. 

In this figure the parabola corresponds to the locus of points satisfying 


<tj ♦ 4^ “ 0 


(15) 
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Within the lower (shaded) region the ACF’s are more complicated. The 
parameters will plot in this region whenever the roots of the equation 

4>(B) - 1 - ♦jB - $ 2 B 2 - 0 (16) 

are complex. When this occurs the process will display a pseudo-periodic 
behavior. 

Not all values of <J>j and are likely to be encountered in this 
study even within the admissible region (Figure 4)* It is probably 
reasonable to assume that as constants of proportionality representing 
the effects of slope wash and creep, the values of a and b ought to 
express the proportional effects of these processes and thus fall in the 
range 0% to 100% (i.e. 0.0 to 1.0). If so then only the region shaded in 
figure 6 should be observed to occur in slope series. An even greater 
constraint will apply if we assume that a and b together should not provide 
more than 100% of the proportional effect, that is: 

a + b < 1 (17) 

It is quite conceivable that they could provide less than 100% of the 
proportional effect since we allow other agents (tectonics, mass wasting, 
climate change, etc.) as represented by A(X) to also have an influence on 
slope forms. Thus, it would appear that we should only expect parameter 
values within the area shown in figure 7. The extent to which this 
turns out to be the case is thus a measure of the appropriateness and 
correctness of the slope process/form relation stated above. 

Another interesting point is to observe that specific values of 
and $2 on ce estimated can be used to predict values of E(X). The 
quality of the prediction to be expected can be measured and is described 
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by the percent variance (of E(X)> explained by the nodal 

E(X) • E (X-l) + <J> 2 E(X-2) 

the remaining variance ia attributed to the additional error term A(X) . 
Thus, the variance of E(X) depends upon that of A(X) according to the 
equation 


^ — Pj4*^ ” P2^2 


(18) 


or, strictly in terms of 0's 


1 “ 


1 + <p 2 I {(1-4)2) ' v 


(19) 


Thus, exactly as would be expected, the percent variance explained by 
these processes increases as the proportional effect increases and is 
zero when their effect is zero. Contours of percent variance explained 
are shown in figure 8. As can be seen, the proportional effects of the 
processes must be quite large before the percent variance explained gets 
close to 100%. 

The final point to consider is the effect upon estimates of 0^ and 4>2 
of errors in collecting data on slope. If the added error constitutes a 
white noise series, W(X) uncorrelated with the original series (in 
particular, uncorrelated with the A(X)) then the resulting series will be 


Z(X) - E(X) + W(X) 


( 20 ) 


The resulting model will be of the form (Box and Jenkins, 1970, p. 121) 

0(B) V d Z(X) - A(X) + 0(B) V^X) (21) 
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Where $(B) le the same polynomial ae given In equation 16 and V* follow 
the definition of difference operator given earlier. Thie ylelda a mixed 
ARIMA model of order (2,1,3). Thua, It ahould be very obvloua If any 
algnlf leant errors are introduced during the data collection proceaa. 

Such a process would have an autocovariance function which would be 
the sum of the autocovariance of the original series plus that of the 
series 


V(X) - V d W(X) 

Y(Z) - Y(E) + Y(V) 


( 22 ) 


and y(V) would consist of a value of -0.5 at lag 1 and would be zero 
everywhere else. Thus, it could be expected that whenever significant 
errors are introduced in data collection the ACF will show the error. 
The value of lag one will be small (less than 0.5) and the ACF will 
otherwise show a normal decay from a value 0.5 greater than that. 
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III. METHODS 


The phyaiographic provincaa choaan for this study ara Indlcatad In 
figura 1 and Tabla I. Two provlncas which could hava baan Included but 
which wara not dua to the limitation to a total of eight ware tha 
Coastal Plain and the Central Lowland. These seamed most reasonable to 
avoid in this first study because it seemed possible that they tend to 
be more strongly affected by depositional processes. Thus, there was 
some uncertainty whether the ARIMA (2,1,0) creep/slope wash model would 
apply. In the case of the Coastal Plain there is a great deal of flood 
plain formation in addition to the marine deposition (and erosion) 
processes. The Central Lowland has so recently been glaciated that the 
deposition of thick sequences of till can be expected to exert a signifi- 
cant control upon the form of the land. The effects of normal erosional 
processes may not be a significant factor in shaping the land yet. 

Choice of specific sampling sites within the chosen physiographic 
provinces was limited by four constraints. Within each region three 
sites were desired. It was assumed that the sites should represent, to 
the extent practical, the diversity of physiography to be found within 
the province. An important guide used to help in this choice was the 
published list of 100 diverse physiographic sites throughout the U.S. 
Whenever possible, sites were chosen from this list (Upton, 1955). 

The scale at which an area is mapped or at which the map is published, 
will Influence the accuracy (and precision) of the data obtained from 
that map. In the case of topographic information the horizontal accuracy 
is directly controlled. The vertical accuracy is also influenced because 
the density of contours is limited, hence so is the contour interval. 

This in turn will control the vertical accuracy (U.S.C.S., 1970) . For 
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this study, only naps of one scale (1:24000) were used. Thus, only 
jhere standard 7 ) 1 ' quadrangle maps were available could sites be selected 
for analysis. 

A further guide to selection of sites Is tsken from the detailed 
breakdown of physiographic provinces into sections. These are contiguous 
regions of similar physiography to some extent distinctive from other 
sections in the province. Sections Into which the eight provinces of 
interest are divided are listed In table II and are shown by light lines 
on figure 1. Wherever possible the three sites within a province were 
taken from three different sections. Four of these provinces have only 
two sections. 

a. Piedmont Province 

b. Blue Ridge Province 

c. Ozark Plateaus Province 

d. Ouachita Province 

If more sites were still available than could be used the selection 
was made randomly. There were also a number of sites that could not be 
selected using the methods outlined above. These were chosen more or 
less arbitrarily. Three constraints wereapplied to this selection. 

The 7V map had to be readily available, the site could not be urbanized 
and the site could not be predominantly in an alluvial plain or other 
depositional setting. These additional sites are listed in table III. 

Three-dimensional plots of the chosen sites showing an area of three 
by three quadrangles surrounding the sampling point were prepared in 
order to examine the areas. This allowed a simple means to familiarize 
oneself with the physiography and to thus ensure that the chosen quadrangle 
was within the province expected and that the region was reasonably 
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Tab It II. Sections into which the eight physiographic 
provincee studied ara divided. 


Physiographic province 


Section 


1. interior Low Plateaus 

2. New England province 

3. Piedmont province 

4. Blue Ridge province 

5. Ozark Plateaus 

6. Valley and Ridge province 

7. Appalachian Plateaus 

8. Ouachita province 


a Highland rim section 
b Lexington Plain 
c Nashville Basin 
d (possible western section) 

a Seaboard Lowland section 
b New England Upland section 
c White Mountain section 
d Green Mountain section 
e Taconlc section 

a Piedmont Upland 
b Piedmont Lowlands 

a Northern section 
b Southern section 

a Sprlngf ield-Salem plateau 
b Boston “Mountains" 

a Tennessee section 
b Middle section 
c Hudson Valley 

a Mohawk section 
b Catsklll section 
c Southern New York section 
d Allegheny Mountain section 
e Kanawha section 
f Cumberland Plateau section 
g Camberland Mountain section 

a Arkansas Valley 
b Ouachita Mountains 


Tabic III. Additional aaaplc aitaa not choacn from 
Upton (1955). 


Quadrangle 

State 

Physiographic province 

Alexandria 

PA 

Valley and Ridge 

Fidelity 

MO 

Ozark Plateaus 

Horseshoe Mtn. 

AK 

Ouachita Mountains 

Lavac « 

AK 

Ouachita Mountains 

Mena 

AK 

Ouachita Mountains 

Saugerties 

NY 

Valley and Ridge 

Sherando 

VA 

Blue Ridge 
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representative of the relief of the chosen physiographic section. 

Earlier investigations had used four traverses oriented at 45° to 
one another along the principal compass directions. This yields replicates 
of the measure of physiographic complexity within each quadrangle which in 
turn allowa us to test whether it varies significantly from site to site 
within a province. These orientations should also be adequate to detect 
the greatest diversity of topographic structure. Thus, this same sampling 
plan was continued in this study. 

These earlier studies had consisted of traverses of points spaced 
2 mm (» 48 m) apart for a distance of 50 points (• 2400 m) . This did not 
seem adequate to ensure covering a sufficiently long part of the area to 
obtain a significant portion of the physiographic diversity. Thus, for 
this study a traverse three times as long (7200 a) was obtained in each 
of the four directions. 

The traverses were sampled at 1 m (■ 24 m) because it was felt 
the earlier samples may have failed to pick up all of the available 
topographic variation. It was assumed that this was the practical limit 
of resolution and the horizontal map accuracy hardly warrants any tighter 
sampling. 

Elevations were estimated to the nearest foot* and recorded b/ 
hand on individual data sheets. Also recorded was the operator's name, 
beginning and ending times, name and location of quadrangle, contour 
interval and any problems that arose including errors in the maps them- 
selves. The traverse positions were determined through the use of a 
transparent overlay (Figure 9) which was centered on the quadrangle. 

Both were taped down on a light table to facilitate the work. All raw 
data are on file with the author. 


* Contours are given in feet on the maps and so metric units were not recorded 


Figure 9. Sampling grid used in collect in traverse 

data. This is a reproduction at about one- 
third original scale. 
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Once recorded the date were keyed into the HP9845B minicomputer end 

stored on floppy disk. Initial accuracy checks included proof-reading 

the data followed by computation of basic statistics including Fishers 

2 

(1953) gj and gj and the X measure of goodness of fit to the normal 
distribution which in most cases appeared tc be the most reasonable 
assumption of the underlying frequency distribution. A typical histogram 
is shown in figure 10, these were produced for each traverse and examined 
carefully for outliers that could indicate erroneous data. A number of 
points were corrected in this way. Summary statistics of raw traverses 
of elevation can be found in the appendix. Data were easily corrected 
using the simple programs available on the HP. In some cases it was 
necessary to refer to the original maps to estimate the correct data. 

In other cases referring to the raw data sheets cleared up the problem. 

Not uncommonly the error was simply a typographical one in which, for 
example, 4160 was keyed in instead of 2160. 

Following this initial data checking the data were transferred by 
phone line to the KSU Burroughs computer system. There the ACF's and 
P ACF's were computed for the original series and the first and second 
differences. In addition, a line-printer generated plot of the data was 
made. This was examined for any obvious points which were erroneous 
values. A significant number of mistakes were caught using this method. 

With the ACF's and PACF's available the next step was to identify 
the model. In this case it was already expected that the model would be 
an ARIMA (2,1,0). Thus, the "identification" step was also a check of 
the reasonableness of the initial assumption. The need for a first 
difference in the model was checked by the decay of the ACF of the original 
series. A slow decay is a good indication that a difference Is needed. 



NORRIS TENNESSEE - TRRVERSE b 



ELEVRTION (feet) 
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The ACF of the first difference was then checked to see that it followed 
the patterns shown in figure 2, parts A or B. The PACF of the first 
difference was checked to see if it followed the desired pattern (figure 3). 
Only two values should be significant, those at lags 1 and 2. The first 
should be positive, the second negative. Any departures from the expected 
pattern were noted and, if appropriate, alternative nodels were identified. 

In every case the ARIMA (2,1,0) model seemed to be the appropriate 
one or at the very least (in a very few cases) was a sufficiently reason- 
able null hypothesis. For each series preliminary estimates of and 

A A 

$2 were obtained from the listed values of and P 2 (see Table IV). 

Whenever the preliminary estimates did not correspond to the expected 

A 

model - this was especially evident if ^ was positive but also could be 
detected if ^ was small (less than 0.5) - the data were checked for 
possible errors. A number of corrections were made at this stage. 

Using these preliminary estimates the actual values of the coefficients 
were obtained (Table V) using the least squares method mentioned pre- 
viously. Confidence limits on those values were also obtained. In 
most cases the preliminary and final estimates corresponded quite closely 
(Figure 11). Where these values did not correspond the original data 
were checked to determine if this were due to an error. This was found 
to be the case in most instances. However, some unexplained large 
discrepancies still remain. Thus, one should not rely too heavily upon 
the preliminary estimates when accurate estimates of and 4*2 are 
essential. 

The reasonableness of the computed model was tested using several 
diagnostic tools. First the residuals, R, were computed by subtracting 
values predicted using this model (with the coefficients 4 >j and $2 


Table iv. Initial estimates of autocorralatlon and preliminary 

coefficient aatlmataa derived from them using equation 13. 


Province 

Quadrangle 

State 

1 

VPi 

1 

V 1 Prelim, 
♦l 

V* Prelim. 
♦ 2 

Interior 

Mammoth 

KY 

a 

.7597 

.5481 

.8119 

>.0687 

Low 

Cave 


b 

.7073 

.4343 

.8007 

-.1320 




c 

.6965 

.4370 

.7616 

-.0934 




d 

.6937 

.3952 

.8087 

-.1658 


Hillsboro 

KY 

a 

.6735 

.3460 

.8061 

-.1969 




b 

.5582 

.1809 

.6642 

-.1898 




c 

.6623 

.3123 

.8114 

-.2251 




d 

.5016 

.1498 

.5667 

-.1333 


Rover 

TN 

a 

.6543 

.3700 

.7208 

-.1016 




b 

.7535 

.5620 

.7635 

-.0133 




c 

.7208 

.5477 

.6786 

.0586 




d 

.7410 

.5679 

.7101 

.0417 


New England 

Ayer 

MA 

a 

.5069 

.2694 

.4984 

.0168 




b 

.6711 

.3609 

.7804 

-.1628 




c 

.8175 

.6998 

.7399 

.0949 




d 

.6772 

.4428 

.6970 

-.0292 


Kingston 

RI 

a 

.6647 

.3978 

.7171 

-.0789 


* 


b 

.4549 

.1794 

.4707 

-.0347 




c 

.6663 

.3349 

.7970 

-.1961 




d 

.5840 

.2384 

.6750 

-.1558 


Brandon 

VT 

a 

.7764 

.6106 

.7611 

.0196 




b 

.6756 

.5740 

.5295 

.2163 




c 

.7075 

.4735 

.7458 

-.0542 




d 

.3504 

.4447 

.6729 

.0522 


Piedmont Warm 

GA 

a 

.6835 

.3384 

.8487 

-.2417 

Springs 


b 

.7249 

.5114 

.7464 

-.0297 


c 

.6341 

.3677 

.6706 

-.0575 



d 

.6689 

.4111 

.7129 

-.0657 

Patterson 

NJ 

a 

.7544 

.5731 

.7474 

.0092 



b 

.7477 

.5046 

.8400 

-.1235 



c 

.6205 

.3077 

.6985 

-.1257 



d 

.7515 

.5397 

.7948 

-.0576 

Washington 

DC 

a 

.5338 

.2092 

.5903 

-.1059 

West 


b 

.6363 

.3440 

.7014 

-.1023 



c 

.6733 

.3579 

.7908 

-.1746 



d 

.6792 

.3196 

.8579 

-.2631 
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v‘p, 

v‘p 2 

V 1 Prelim. 

V 2 Prelim. 

Province 

Quadrangle 

State 

♦» 

♦ 

2 

Blue Ridge 

Mount 

NC 

a 

.8078 

.6789 

.7465 

.0759 

Mitchell 


b 

.8430 

.7278 

.7930 

.0593 




c 

.6686 

.3641 

.7689 

-.1500 




d 

.8226 

.6811 

.8113 

.0137 


Strasburg 

VA 

a 

.7714 

.6113 

.7377 

.0420 



b 

.65 

.33 

.7541 

-.1602 




c 

.73 

.54 

.7189 

.0152 




d 

.81 

.66 

.8088 

.0113 


Sherando 

VA 

a 

.8600 

.77 

.7596 

.1167 




b 

.87 

.76 

.8589 

.0128 




c 

.8386 

.7163 

.8017 

.0440 




d 

.8381 

.7066 

.8274 

.0149 


Ozark 

Plateaus 


I ronton 


Saint Paul 


Fidelity 


MO a 
b 


AR a 
b 


MO a 
b 


Ouachita 

Mtns. 


Horseshoe 

Mtn. 


AR a 
b 


AR a 
b 


AR a 

b 


.7758 

.6144 

.7514 

.0315 

.8096 

.6455 

.8330 

-.0289 

.8063 

.7079 

.6731 

.1651 

.7359 

.5373 

.7427 

-.0093 

.7972 

.6392 

.8018 

-.0171 

.8401 

.7342 

.7589 

.0966 

.6997 

.4434 

.7630 

-.0905 

.7335 

.6157 

.6102 

.1681 

.65 

.19 

.9117 

-.4026 

.58 

.33 

.5856 

-.0096 

.43 

.29 

.3746 

.1289 

.6140 

.2679 

.7092 

-.1626 

.7928 

.6232 

.8042 

-.0144 

.7832 

.6114 

.7873 

-.0052 

.8441 

.6903 

.9093 

-.0772 

.8074 

.6800 

.7422 

.0807 

.6807 

.4062 

.7532 

-.1065 

.6500 

.3700 

.7081 

-.0909 

.5500 

.2400 

.5993 

-.0896 

.64 

.32 

.7371 

-.1518 

.77 

.49 

.9646 

-.1528 

.73 

.66 

.5314 

.2721 

.75 

.63 

.6343 

.1543 

.75 

.53 

.8057 

-.0743 


Lavaca 
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Province 

Quadrangle 

State 


v‘p 2 

V 1 Prelim. 
♦, 

V 2 Prelim. 
♦ 

2 

Valley and 

Norris 

TN 

a 

.7677 

.3667 

.8101 

-.0552 

Ridge 



b 

.7319 

.5092 

.7736 

-.0570 




c 

.8167 

.6251 

.9195 

-.1258 




d 

.7896 

.5610 

.9206 

-.1659 


Alexandria 

PA 

a 

.7481 

.5204 

.8148 

-.0891 




b 

.6833 

.4332 

.7265 

-.0632 




c 

.7807 

.5376 

.9244 

-.1841 




d 

.7071 

.4970 

.7113 

-.0060 


Saugerties 

NY 

a 

.6427 

.2760 

.7928 

-.2335 




b 

.7318 

.4968 

.7928 

-.0834 




c 

.7186 

.3846 

.9144 

-.2725 




d 

.6330 

.3258 

.7121 

-.1250 


Appalachian Ithaca NY a 

Plateaus West b 

c 

d 

Fayetteville WV a 

b 
c 
d 

TN a 
b 
c 
d 


,7323 

.5201 

.7578 

-.0349 

,6572 

.2572 

.8593 

-.3075 

,5511 

.8029 

.1560 

.7169 

7168 

.4288 

.8421 

-.1748 

6652 

.5825 

.4985 

.2560 

8215 

.6916 

.7792 

.0515 

7368 

.6223 

.6088 

.1738 

7106 

.5351 

.7027 

.0241 

8977 

.8287 

.8053 

.1053 

4471 

.2816 

.4014 

.1021 

6981 

.4376 

.7658 

-.0970 

8060 

.6743 

.7493 

.0704 


Whitwell 
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Table V. Final estimates of coefficients and $2 * or 96 traverses. 

Also listed are the upper and lower 95% confidence limits for the 
coefficients. 

01 01 02 02 
Province Quadrangle State Upper Lower Lower 02 Upper 


Piedmont Warm Springs GA a 

b 

c 

d 

Paterson NJ a 

b 

c 

d 

Washington DC a 

West b 

c 
d 


Blue Ridge Mount NC a 

Mitchell b 

c 
d 

Strasburg VA a 

b 

c 

d 

Sherando VA a 

b 

c 

d 


Valley and Norris TN a 

Ridge b 

c 
d 

PA a 
b 
c 
d 


0.97 

0.86 

0.74 

- 0.35 

- 0.24 

- 0.13 

0.86 

0.74 

0.63 

- 0.14 

- 0.03 

0.09 

0.79 

0.67 

0.55 

- 0.17 

- 0.06 

0.06 

0.83 

0.71 

0.60 

- 0.18 

- 0.06 

0.05 

0.93 

0.81 

0.69 

- 0.19 

- 0.07 

0.05 

0.96 

0.84 

0.73 

- 0.24 

- 0.12 

- 0.05 

0.74 

0.63 

0.53 

- 0.18 

- 0.08 

0.03 

0.92 

0.80 

0.68 

- 0.17 

- 0.05 

0.07 

0.71 

0.59 

0.48 

- 0.22 

- o.n 

0.01 

0.82 

0.71 

0.60 

- 0.22 

- o.n 

0.01 

0.91 

0.79 

0.68 

- 0.29 

- 0.17 

- 0.06 

0.97 

0.86 

0.75 

- 0.38 

- 0.26 

- 0.15 

0.86 

0.75 

0.63 

- 0.05 

0.08 

0.19 

0.91 

0.80 

0.68 

- 0.05 

0.06 

0.18 

0.89 

0.77 

0.66 

- 0.26 

- 0.15 

- 0.03 

0.93 

0.81 

0.69 

- 0.10 

0.02 

0.13 

0.86 

0.75 

0.63 

- 0.07 

0.04 

0.16 

0.87 

0.75 

0.64 

- 0.27 

- 0.15 

- 0.04 

0.84 

0.72 

0.61 

- 0.10 

0.01 

0.13 

0.93 

0.81 

0.70 

- 0. 10 

0.02 

0.14 

0.89 

0.77 

0.66 

- 0.01 

0.11 

0.22 

0.98 

0.86 

0.74 

-0.10 

0.02 

0.14 

0.92 

0.81 

0.69 

- 0.07 

0.05 

0.16 

0.94 

0.83 

0.71 

- 0.10 

0.02 

0.13 

0.93 

0.81 

0.70 

- 0.17 

- 0.06 

0.06 

0.91 

0.79 

0.68 

- 0.18 

- 0.06 

0.05 

1.03 

0.92 

0.80 

- 0.24 

- 0.12 

- 0.01 

1.03 

0.92 

0.80 

- 0.28 

- 0.16 

- 0.05 

0.93 

0.82 

0.70 

- 0.20 

- 0.09 

0.03 

0.84 

0.73 

0.61 

- 0.17 

- 0.06 

0.06 

1.05 

0.93 

0.82 

- 0.30 

- 0.19 

- 0.08 

0.83 

0.71 

0.59 

- 0.12 

-0.01 

0.11 


Alexandria 
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Province 

Quadrangle 

State 

Upper 

♦l 

Lower 

Lower 

^2 

Upper 

Valley and 

Saugertlea 

NY a 

0.91 

0.79 

0.68 

-0.35 

-0.23 

-0.12 

Ridge 


b 

0.91 

0.79 

0.68 

-0.20 

-0.08 

0.03 



c 

1.05 

0.93 

0.82 

-0 39 

-0.28 

-0.17 



d 

0.81 

0.70 

0.59 

-0.23 

-0.12 

-0.01 


Appalachian 

Ithaca West 

NY 

a 

0.91 

0.79 

0.68 

-0.12 

-0.04 

0.11 

Plateaus 



b 

1.00 

0.88 

0.77 

-0.41 

-0.30 

-0.18 




c 

0.84 

0.73 

0.62 

-0.43 

-0.32 

-0.21 




d 

1.05 

0.94 

0.82 

-0.28 

-0.17 

-0.06 


Fayetteville 

WV 

a 

0.61 

0.50 

0.39 

0.14 

0.25 

0.36 




b 

0.90 

0.78 

0.66 

-0.06 

0.05 

0.17 




c 

0.72 

0.61 

0.49 

0.06 

0.17 

0.29 




d 

0.78 

0.67 

0.55 

-0.05 

0.06 

0.18 


Whitwell 

TN 

a 

0.89 

0.77 

0.66 

0.02 

0.14 

0.25 




b 

0.51 

0.40 

0.28 

-0.01 

0.11 

0.22 




c 

0.89 

0.78 

0.66 

-0.21 

-0.09 

0.02 




d 

0.87 

0.75 

0.63 

-0.04 

0.07 

0.39 


Interior Mammoth Cave 

KY 

a 

0.98 

0.86 

0.75 

-0.22 

-0.11 

0.01 

Low 


b 

0.92 

0.80 

0.68 

-0.25 

-0.13 

-0.02 



c 

0.89 

0.78 

0.66 

-0.21 

-0.09 

0.02 



d 

0.95 

0.84 

0.72 


-0.19 

-0.07 

Hillsboro 

KY 

a 

0.91 

0.80 

0.68 


-0.19 

-0.08 



b 

0.77 

0.65 

0.54 

-0.30 

-0.18 

-0.07 



c 

0.92 

0.81 


-0.34 

-0.22 

-0.11 



d 

0.69 

0.57 

0.46 

-0.25 

-0.14 

-0.02 

Rover 

TN 

a 

0.84 

0.73 

0.61 

-0.23 

-0.11 

0.01 



b 

0.88 

0.76 

0.65 

-0.12 

-0.01 

0.11 



c 

0.80 

0.68 

0.56 

-0.06 

0.06 

0.18 



d 

0.83 

0.71 

0.59 

-0.07 

0.04 

0.16 


New England Ayer 

MA 

a 

0.61 

0.50 

0.38 

-0.10 

0.02 

0.13 



b 

0.90 

0.78 

0.67 

-0.28 

-0.16 

-0.05 



c 

0.88 

0.77 

0.65 

-0.03 

0.09 

0.21 



d 

0.81 

0.70 

0.58 

-0.14 



Kingston 

RI 

a 

0.85 

0.73 

0.61 

-0.19 

-0.07 

0.05 



b ' 

0.56 

0.45 

0.33 

-0.13 

-0.02 

0.08 



c 

0.69 

0.57 

0.46 

-0.17 


0.07 



d 

0.79 

0.68 

0.56 

-0.27 

-0.15 

-0.04 

Brandcn 

VT 

a 


0.79 

0.67 

-0.10 

0.02 

0.14 



b 

0.66 

0.54 

0.43 

0.12 

0.23 

0.35 



c 

0.86 

0.75 

0.64 

-0.16 


0.06 



d 

0.51 

0.39 

0.27 


-0.08 

0.04 
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r 


Province 


Otark 

Plataaua 


Ouachita 

Mountains 


Quadrangle 

State 


♦l 

*1 

Lower 

92 

Lower 

♦ 2 

92 

Upper 

I ronton 

MO 

a 

0.87 

0.76 

0.64 

-0.08 

0.04 

0.15 



b 

0.95 

0.83 

0.72 

-0.14 

-0.03 

0.09 



c 

0.79 

0.67 

0.56 

0.05 

0.16 

0.28 



d 

0.86 

0.74 

0.63 

-0.12 

-0.001 

0.12 

Saint Paul 

AR 

a 

0.93 

0.82 

0.70 

-0.13 

-0.02 

0.09 



b 

0.88 

0.76 

0.65 

-0.02 

0.01 

0.21 



c 

0.88 

0.76 

0.65 

-0.20 

-0.09 

0.03 



d 

0.79 

0.68 

0.56 

0.01 

0.12 

0.23 

Fidelity 

MO 

a 

1.01 

0.90 

0.79 

-0.50 

-0.39 

-0.28 



b 

0.78 

0.65 

0.53 

-0.16 

-0.04 

0.08 



c 

0.50 

0.38 

0.27 

0.02 

0.13 

0.25 



d 

0.84 

0.72 

0.61 

-0.29 

-0.17 

-0.06 


Horseshoe 

AR 

a 

0.92 

0.80 

0.69 

-0.12 

-0.01 

0.11 

Mtn. 


b 

0.93 

0.82 

0.70 

-0.14 

-0.02 

0.10 



c 

1.01 

0.90 

0.78 

-0.18 

-0.06 

0.05 



d 

0.86 

0.74 

0.63 

-0.03 

0.08 

0.20 

Mena 

AR 

a 

0.87 

0.76 

0.64 

-0.22 

-0.10 

0.02 



b 

0.82 

0.71 

0.59 

-0.19 

-0.07 

0.04 



- 

0.48 

0.59 

0.70 

-0.19 

-0.07 

0.04 



d 

0.87 

0.76 

0.64 

-0.27 

-0.16 

-0.49 

Lavaca 

AR 

a 

1.07 

0.96 

0.84 

-0.36 

-0.24 

-0.13 



b 

0.64 

0.52 

0.41 

0.18 

0.29 

0.40 



c 

0.72 

0.60 

0.49 

0.08 

0.19 

0.30 



d 

0.92 

0.80 

0.69 

-0.19 

-0.07 

0.04 
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listed in Table V) from the values actually observed: 

E*(X) - ifrjE^X-l) - 4> 2 E 1 (X-2) » R(X) (23) 

The residuals should be white noise, i.e. have no remaining structure. 

This was tested by computing the ACF’s and FACF's as previously. The 

ACF and PACF of the original (residual) series should show no significant 

values. The ACF of the first difference is expected to have a value of 

-0.5 at lag 1 and be not significantly different from zero otherwise. 

The PACF of the first difference is expected to decay to zero with the 

value at each lag k equal to -l/(k+l). The ACF of the second difference 

should have the values -0.67 at lag 1 and 0.17 at lag 2 and otherwise be 

not significantly different from zero. The PACF of the second difference 

should decay to zero following the function -2/(k+2). Each of these 

checks wore performed and significant departures recorded. 

An additional diagnostic is that the mean of the residual series 

should not be significantly different from zero. This was tested by a 

t statistic (Table VI). The "portmanteau" test of the ACF of the original 

residual series checks for systematic tendencies to deviate from zero 

which is more sensitive than simply checking the confidence limits. This 

2 

test yields a statistic distributed as X under the null hypothesis that 
there are no significant systematic tendencies to deviate from zero. The 
results of this test are also reported in table VI. 

The fit of the model was also checked by examining a plot of the 
residuals (versus observation number). Suspicious patterns can be an 
indication of lack of fit. Such checks were also done, a typical plot 
of residuals is given in figure 12. 

Finally, because it was expected that the coefficients should fall 
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Table VI . Results of first five tests of residuals from the ARIMA 
(2,1,0) model. Tests are explained in the text. An 
asterisk indicates the value exceeds the 95Z confidence 
limit. 


Province 

Quadrangle 

State 

x/st. er. 

Y 2 

x 46 

v'pi 

V 2 Pl 

v 2 p 2 

Interior 

Mammoth 

KY 

a 

.07 

57.3 

-.56 

-.71 

.24 

Low 

Cave 



.26 

43.8 

-.49 

-.67 

.22 




c 

1.10 

35.8 

-.50 

-.65 

.11 




d 

.01 

52.9 

-.52 

-.69 

.23 


Hillsboro 

KY 

e 

.02 

58.1 

-.57 

-.73 

.28 




b 

.40 

72.8* 

-.53 

-.71 

.27 




c 

.16 

49.3 

-.52 

-.68 

.19 




d 

.35 

72.3* 

-.57 

-.75 

.38* 


Rover 

TN 

a 

.13 

40.4 

-.55 

-.73 

.35* 




b 

.76 

40.6 

-.52 

-.67 

.15 




c 

.44 

39.1 

-.50 

-.67 

.18 




d 

.01 

55.1 

-.50 

-.65 

.13 

Mew England 

Ayer 

MA 

a 

.11 

68.5* 

-.50 

-.65 

.07 




b 

.25 

67.3* 

-.51 

-.69 

.24 




c 

.27 

67.5* 

-.52 

-.67 

.15 




d 

.05 

44.8 

-.49 

-.65 

.10 


Kingston 

RI 

a 

1.57 

38.8 

-.53 

-.71 

.32 




b 

.08 

67.5* 

-.49 

-.67 

.24 




c 

.25 

45.6 

-.54 

-.73 

.35* 




d 

.60 

84.4* 

-.51 

-.68 

.20 


Brandon 

Vi 

a 

1.63 

64.5* 

-.53 

-.69 

.18 




b 

1.78 

55.8 

-.47 

-.66 

.20 




c 

.10 

28.7 

-.50 

-.68 

.21 




d 

1.18 

29.5 

-.46 

-.61 

.03 

Piedmont 

Warm 

GA 

a 

1.14 

44.2 

-.54 

-.69 

.20 


Springs 


b 

.14 

42.0 

-.50 

-.65 

.13 




c 

.09 

44.1 

-.54 

-.70 

.24 




d 

.78 

46.0 

-.53 

-.69 

.19 


Paterson 

NJ 

a 

.55 

66.8* 

-.53 

-.70 

.73* 




b 

.86 

68.9* 

-.46 

-.62 

.02 




c 

1.56 

33.8 

-.49 

-.66 

.15 




d 

1.03 

41.5 

-.48 

-.64 

.09 


■ ■— 

Province 


Piedmont 


Quadrangle 


Washington 

West 


State 

x/st. er. 

X 

M 

-r* 

ON 


V 2 Pl 

v 2 p 2 

DC a 

.02 

60.1 

-.46 

-.63 

.09 

b 

.19 

44.1 

-.47 

-.64 

-.13' 

c 

.83 

50.3 

-.51 

-.67 

.19 

d 

.87 

79.1* 

-.47 

-.64 

.12 




Blue Ridge 




I 


Mount 

NC 

a 

.10 

23.2 

-.52 

-.69 

.23 

Mitchell 


b 

.98 

39.9 

-.49 

-.66 

.17 



c 

1.08 

45.9 

-.49 

-.65 

.11 



d 

.35 

50.3 

-.52 

-.69 

.20 

Strasburg 

VA 

a 

1.07 

39.3 

-.51 

-.67 

.16 



b 

.90 

63.3* 

-.51 

-.67 

.14 



c 

.41 

48.2 

-.50 

-.68 

.22 



d 

.33 

31.5 

-.48 

-.64 

.13 

Sherando 

VA 

a 

.73 

43.3 

-.51 

-.67 

.16 



b 

.72 

54.2 

-.49 

-.63 

.06 



c 

1.30 

45.8 

-.49 

-.65 

.16 



d 

.40 

76.3* 

-.49 

-.66 

.13 


Ozark 

Plateaus 


Ironton 

MO 

a 

1.14 

58.8 

-.51 

-.66 

.15 




b 

.56 

46.4 

-.54 

-.70 

.24 




c 

.13 

47.5 

-.50 

-.67 

.18 




d 

.26 

58.6 

-.45 

-.64 

.16 

l 

Saint Paul 

AR 

a 

.05 

35.6 

-.48 

-.64 

.08 

j 



b 

.41 

47.7 

-.53 

-.70 

.26 




c 

.79 

80.3* 

-.54 

-.71 

.25 

j 



d 

.08 

60.7 

-.54 

-.72 

.29 

■] 

Fidelity 

MO 

a 

.50 

58.8 

-.60 

-.75 

.36* 

.! 



b 

1.80 

36.0 

-.52 

-.72 

.33 




c 

1.32 

100.7* 

-.49 

-.64 

.08 

j 



d 

.16 

31.7 

-.51 

-.67 

.18 



Ouachita Horseshoe 

Mtns. Mtn. 


AR 

a 

.68 

69.3* 

-.52 

-.69 

.26 


b 

.56 

82.4* 

-.49 

-.64 

.08 


c 

1.16 

68.4* 

-.59 

-.72 

.24 


d 

.47 

61.2 

-.5? 

-.73 

.29 

AR 

a 

1.46 

41.3 

-.49 

-.65 

.10 


b 

1.53 

44.5 

-.49 

-.66 

.15 


c 

.53 

30.3 

-.50 

-.67 

.13 


d 

.65 

25.4 

-.47 

-.64 

.13 


Mena 
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Province 

Quadrangle 

State 

x/at. er. 

2 

X 46 

V Pi 


** 

Ouachita 

Lavaca 

AR a 

.49 

86.4* 

-.62 

-.78 

.43* 

Mtns. 


b 

1.02 

39.6 

-.51 

-.67 

.18 



c 

1.02 

37.3 

-.45 

-.64 

.14 



d 

1.06 

81.0* 

-.50 

-.66 

.17 


Valley and Norris 

TN 

a 

.08 

34.3 

-.51 

-.65 

.14 

Ridge 


b 

.79 

51.4 

-.56 

-.73 

.32 



c 

.25 

58.0 

-.46 

-.60 

-.01* 



d 

.08 

54.6 

-.49 

-.65 

.10 

Alexandria 

PA 

a 

.87 

53.1 

-.50 

-.67 

.18 



b 

.61 

23.5 

-.51 

-.66 

.18 



c 

.60 

41.8 

-.47 

-.65 

.14 



d 

.10 

70.3* 

-.42 

-.57 

-.06* 

Saugerties 

NY 

a 

.14 

42.9 

-.44 

-.60 

.01 



b 

.07 

41.3 

-.47 

-.63 

.08 



c 

1.50 

83.3* 

-.47 

-.63 

.04 



d 

.47 

45.8 

-.48 

-.68 

.09 


Appalachian 

Ithaca West 

NY 

a 

2.77 

61.4 

-.41 

-.58 

0.0* 

Plateaus 



b 

2.78 

35.2 

-.51 

-.69 

.24 




c 

1.54 

62.6 

-.55 

-.71 

.23 




d 

3.53 

38.2 

-.50 

-.69 

.26 


Fayetteville 

WV 

a 

.31 

73.6* 

-.47 

-.64 

.15 




b 

.00 

75.1* 

-.47 

-.65 

.17 




c 

.21 

63.0* 

-.53 

-.72 

.32 




d 

.25 

82.4* 

-.46 

-.64 

.12 


Whitwell 

TN 

a 

.06 

71.7* 

-.45 

-.64 

.15 




b 

.45 

54.8 

-.48 

-.64 

.08 




c 

1.50 

57.6 

-.47 

-.66 

.19 




d 

.40 

40.7 

-.50 

-.66 

.14 


WRSHINGTON WEST 




Figure 12. Typical plot of residuals from the best fit ARINA (2,1,0) 
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within certain limits, as discussed previously, any deviations from 
expectation were taken as an Indication of possible errors In the data. 
This was found to be the case for 18 traverses. Detection of the actual 
error was more difficult in these cases since the obvious ones had much 
earlier been detected. A typical error is illustrated in figure 13, 

Most of those errors were found to be typographical in nature. A special 
program was set up to examine traverses to detect errors such as that in 
figure 13. This program is listed in the Appendix. Where errors such 
as that in figure 13 could not be found, suspicious traverses were checked 
point by point against the original data sheets. If an error still was 
not found the original maps were again referred to. Other possible 
sources of error will be mentioned in Section VII. 

Once the fits and estimates were accepted as reasonable the values 
of and $2 were plotted for each quadrangle (Appendix) . This gives an 
idea of the uniformity of values within the quadrangle and is another 
means of detecting inaccurate estimates. When one value deviated greatly 
from the other three, even though it still fell within the acceptable 
area, the original traverses were again checked for possible erroneous 
values. On the triangular plots of versus $2 are also listed a number 
of summary statistics from the analyses for that quadrangle as well as a 
summary of the climatic data for the nearest station and information 
about the quadrangle map itself . 

Once the data had been thus thoroughly checked and subject to 
scrutiny they were subject to analysis to test the concepts presented in 
the Introduction. A separate Analysis of Variance was done for each 
parameter ($j and (J^). The design was a nested one with quads nested 
within province. Traverses, with four levels, were considered crossed 
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with province in order to consider the possibility thst there is a 
systematic tendency for say slopes on north-south oriented traverses to 
contain a different process mix than other slopes. Following the analysis 
of variance a multiple comparison test was performed to determine which 
levels of the significant factors could be considered significantly 


distinct 
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IV. RESULTS 

This chapter will be divided Into two parte in the firat the appro- 
prlateneaa of the original aeeumptlona about proceea representation by 
meant of an ARIMA (2,1,0) model will be evaluated through a number of 
tests. Results presented in that part appear to fully substantiate the 
chosen model. Part two therefore uses the estimated paremeters of that 
model to examine the question of significant variation in process rates 
from quadrangle to quadrangle within a physiographic province and from 
province to province. In this part the results of the analysis of 
variance and multiple comparison tests are presented. 

Fortunately, there are a large number of tests available to check the 
appropriateness of the hypothesized process form association. They may 
be conveniently divided into two groups. The first addresses the ques- 
tion: How likely is the ARIMA (2,1,0) model to be correct compared to 

other alternatives? The second group examines the goodness of fit of the 
actual estimated parameters since they are further constrained by geomorphic 
considerations. 

Comparison of the ARIMA (2,1,0) to other alternatives can be accom- 
plished satisfactorily by examining the ACF's and PACF's of the elevation 
traverses and their first differences (i.e. slopes). An alternative model 
would still be a member of the ARIMA (p,d,q) family. Thus, they must be 
recognized by showing a different value of p, d or q or some combination 
of these three. Each of these can be detected using the ACF's and PACF's. 

It has been assumed that essentially all traverses under study will 
show the same model (although the parameter values may differ). An 
alternative that ought to be considered is that there is no systematic 
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structure to elevation traverses; any observed structure nay simply be 
coincidental; each area could show a different model and the model (or 
models) may carry absolutely no implication about the processes. It 
should also be remembered that the processes under consideration do not 
comprise the entire set of geomorphlc processes. Certainly mass wasting 
influences landforms, especially in areas of high relief. Thus, some 
departures from the model can be expected, if such areas ere included in 
the study. 

Of the three parameters to be considered, q, can be disposed of most 
readily. If q is non-zero the PACF will show a sinusoidal pattern or a 
damped exponential form or some mixture of the two. This would be the 
case if q was 1 or 2; more complicated models (q > 3) clear ly do not 
app’v to these data. Whether the q 0 pattern shows up in the PACF of 
the original series or in the dth difference series (slopes, changes in 
slope, etc.) depends upon the correct value of d in the model. As will be 
shown, it is highly likely that d is 1 (as also predicted by theory). 

Thus, we can expect that any q 4 0 pattern will show up in the PACF of 
the first differences (i.e. the slope series). These were examined with 
especial interest in the occurrence of exponential or damped sinusoidal 
patterns (Table VII). As can be seen from inspection of the table, such 
patterns are "found" in only eight of the 96 traverses. In each case 
these "patterns" are quite problematical; whether they actually exist is 
open to question. A conservative approach was taken in that anything 
that was anywhere close to the form of interest was accepted as being 
real. Note that none of these so-called patterns actually extend beyond 
the significance limits. The most convincing example is shown in figure 14. 
It is also important to note that comparable patterns were found in the 


Table VIL Occurrence of significant lags and suspected patterns in the PACF of the slopes. In the first coltmm for 

each traverse the numbers indicate which lags (out of the first ten) have significant values. In the pattern 
column Y indicates a simple pattern may exist; see footnote for a description, N indicates there is no simple 
pattern. 
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PACF's of at least five other slope series but disappeared after known 
data errors were corrected. Thus, it is quite possible that these patterns 
may simply be an indication that further (necessarily minor) errors 
remain in these series. At any rate, the information displayed in 
table VII is (in at least 92% of the cases) precisely what would be 
expected if the ARIMA (2,1,0) model is correct. It is the authors opinion 
that no reasonable interpretation of this information would allow general 
rejection of the ARIMA (2,1,0) model. Furthermore, it would appear on 
the basis of this information that the same single model applies nearly 
universally to all of the study areas. 

Determination of the correct value of d is both a critical test of 
the model and conveys significant geomorphic meaning. To claim that d * 1 
implies that the correct phenomena to study is the behavior of, and influ- 
ences upon, slopes rather t an, for example, raw elevations themselves. 

It also implies that elevations are a non-stat ionary phenomenon which has 
also been claimed by Mandelbrot (1977). Examination of d will be divided 
into four parts. First, for each series we will consider if any differencing 

is needed at all (Ho: d * 0; Ha: d > 0) . It will be shown that the null 

hypothesis can be rejected universally in favor of the alternative; some 
kind of differencing should be done. Next, we will consider if perhaps a 
higher order difference is not preferable (Ho: d * 1; Ha: d > 1). It 
will be shown that in a large number of cases this would definitely be 
overdoing it. That is, d is nt most equal to one. Third, we will test 
whether d * 1 is sufficient . The idea is that even if d * 2 could be 

accepted in some cases, it is also the case that d ■ 1 would yield an 

equally acceptable (or more acceptable) model but would be more parsimonious 


and thus preferable. Finally, we will consider the question of the 


I 
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universality of the d « 1 part of the model. 

Three interrelated bits of data are available that show that uae of 
the original data (as opposed to some appropriate difference) is unaccept- 
able. The very high values of (none less than 0.97) in the ACF of the 
original series (Table VIII) imply that without an explicit difference 
(i.e. with d » 0) the series would universally require a model which had 
a ^ parameter nearly or exactly equal to 1.0. This is essentially equal 
to a model with d * 1, it implies a non-s tat ionary model and would make 
estimates of other parameters unreliable. Another strong indication of 
non-stationarity is provided by the very slow rate at which the ACF decays 
from this large initial value. Most of the series show ACF's that remain 
positive throughout the entire range of lags computed (either 25 or 48). 

Of those that do reach aero, the quickest "decay" is in 23 lags (Table VIII) 
although the decay pattern remains strong throughout and no "noisy" random 
fluctuations such as would be expected in a stationary series are observed. 
In general, a stationary model would be expected to display such "noisy" 
behavior by the 10th or, at the very most, by the 20th lag. Thus, the 

rate of decay of the series give a universal indication that differencing 

2 

is needed. Finally, the extremely large values of x obtained from the 
portmanteau test of the ACF of the original series (all exceed 1000 although 
the '■'■% cut-off is 65.17) provide strong indication that the original 
series are not remotely near being white noise (p * d * q ■ 0) and that 
they are almost certainly non-stat ionary . Thus, the combination of the 
three lines of evidence as well as geomorphic reasoning provide a clear 
indication that the original series (elevations) are not appropriate for 
analysis and that some value of d greater than zero will be required to 


produce such a series. 


Table VIII. Data from autocorrelation functions of original series and 
first and second differences relevant to the Identification 
of the appropriate value of d in the ARIMA (p,d,q) model. 


► 


Column Headings 


4 Autocorrelation of original series at lag one. 

5 Lag at which first negative value is found. An N means 
that no negative values were computed. 

6 Number of lags computed. This value is also the number 
of degrees of freedom in the portmanteau tests. 

7 Value of in portmanteau test of ACF of the original 
series. 

2 

8 Value of x in portmanteau test of the first difference. 

2 

9 Ratio of x from original series to that of first difference. 
Less than a sixfold reduction is indicated by an asterisk (*) . 

2 

10 Value of x in portmanteau test of the second difference. 
Significant values (exceeding 65.17) are indicated by 
asterisks (*) . 

11 Significant negative values at lag one of the second 
difference series. A "Y" means significant, and "N" not 
significant. 
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That some differencing is always required is now quite clear. Is it 
first order or higher? The possibility that a higher order difference 
would be preferable is to be considered next. This possibility is tested 
by comparing the results of the portmanteau tests of first and second order 
differences, by examining the sufficiency of a straight second order 
difference model, and by considering the necessity of a second order 
difference. 

Table VIII shows that a first order difference significantly reduces 

the values of x in the portmanteau tests. Whereas the value is always 

above 10C0 for the undifferenced series it is rarely above 1000 when d ■ 1 

(only 6% of the cases) and is usually in the range 100 to 600. These 

latter values are reasonable if the full model contains parameters 

in the range 0.5 to 1.0 as anticipated in the theoretical model (recall 

chapter II, "Theory," page ). The extent that the first order difference 
2 

reduces the x value is also listed in Table VIII. As can be seen, the 
typical reduction is about an order of magnitude, the smallest is by a 
factor of 2.6; only 13% of the reductions are less than a factor of 6. In 
contrast to this, the reduction in going from d * 1 to d * 2 is not nearly 
as drastic. The reduction is rarely as great as a factor of 6. If the 
first difference series were non-stat ionary we could expect at least as 
large a reduction, if not more. As will be shown next in the discussion 
of the p » 2 terms, the ACF of the first difference shows consistent 
evidence of stationarity . Thus, it appears quite definite that all of the 
non-stationarity is removed by a first difference. The only structure 
renaining after the first difference, and producing significant values of 
X in the portmanteau test, is that due to stationary terms (in particular 
$*6, such as and 4>2 as hypothesized in the theoretical model). Thus, 
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the first order difference model Is quite definitely sufficient to ensure 
statior.arity (providing p > 0). 

In contrast it is clear from Table VIII that a simple second order 
model (i.e. p ■ q ■ 0, d ■ 2) is not sufficient . The portmanteau test 
shows that white noise is not acheived in at least 48% of the traverses. 
Thus, if it were true that d ■ 2 it would for at least these traverses be 
necessary to add additional terms (probably 4>'s). Thus, the models would 
be at least as complicated as the one theorized. Additionally, the d ■ 2 
concept would carry with it the implication that models differ from site 
to site and no general geomorphic process/form model can be discerned. 

Not only is the d ■ 2 model not sufficient, it appears that it is not 
necessary in that it provides more model than is appropriate; it overdoes 
the job. This can be seen from the values of in the ACF of the second 
differences. A large proportion (39%) of the series have significant 
negative values at lag 1; a large number of other series show large but 
not significant negative values. Such a pattern is usually taken as prima 
facie evidence that the series has been over-differenced. The resulting 
model would require <J>j < 0 and this essentially is only a correction 
for taking d * 2 rather than d * 1. Althou^ * this particular evidence of 
overdifferencing is not universal, it does occur in at least one traverse 
in all quadrangles but four. Thus, it is again clear that a d ■ 2 idea 
carries with it the implication that no one model is sufficient and that 
these procedures cannot produce geomorphic process/form understanding. 
Although this could at best imply that two models are sufficient, study 
of the remaining four quadrangles provide strong evidence that a great 
variety of different models would be required. Table IX lists the different 
significant lags in these traverses when d * 2. Addition of one parameter 
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certainly cannot produce auch a variety of patterns. Thus, setting d • 2 
could only lead to obfuscation of the relatione which exist between 
adjacent slopes. 

It appears from the above evidence that a first difference is 
certainly necessary and is quite likely sufficient to model the autocor- 
relation of the traverses when combined with some $ terms. It is clear 
that all traverses universally require at least d ■ 1 so that choice of 
this model is satisfying in its general applicability. It remains to be 
seen whether this universality can be maintained in the choice of p. 

Next, we will consider the evidence concerning the appropriateness of 
the p - 2 portion of the model. As shown in chapter II any of four patterns 
in the ACF could be expected if it is true. Geomorphic considerations, 
however, limit the expected parameter values to the area shown in figure 7, 
therefore, only patterns of types A and B (Figure 2) should be found. 

Thus, the ACF should appear to be mixed exponentials or a damped sine 
wave. Presumably near the boundary of the two regions a combination of 
the two could occur. Other patterns would be an indication that the p ■ 2 
idea is not correct. At this point it should be recalled that it is 
reasonable to consider p * 1 to be a legitimate subset of the general 
model. This is so if $2 * 0 which could occur if the landform is not 
significantly shaped by creep. If p ■ 1, the ACF will look like a pure 
exponential decay. The ACF's were examined to determine what sort of 
pattern could be observed. As can be seen in Table X all traverses show 
either an exponential decay* a damped sinusoidal form or a combination of 
the two. Five of the 96 traverses show a very slow decay which could be 
interpreted as a possible need for futher differencing. In no case is 
the decay as slow as was seen in the original series. As will be shown 
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Tabic X. Pattern* observed in the autocorrelation function of the first differences 
of the traverses. 
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later, in all caaaa tha estimated parameter, is below 1.0 and In only 
one of these five cases does the 95% confidence limit about the matliMted 
value even include 1.0. This is strong evidence that further differencing 
is not required. Only one traverse (Brandon Vermont, d) shows a strong 
departure from the basic exponential/sinusoidal pattern. In this traverse 
a mixed model (i.e, both p i 0 and q i 0) is suggested. However, the PACF 
does not support this notion since it does not display the necessary 
exponential/sinusoidal pattern. Thua, it appears that the ARIMA (2,1,0) 
model although a poor fit is apparently the best one available even in 
this case. Interestingly, as Table V shows, this traverse did produce an 
unusually low value for (0.39) whereas almost all other estimates 
exceed 0.50. 

It is concluded on the basis of the patterns observed in the ACF that 
a value of p * 2 is the most reasonable one for the traverses. For some 
it could be found that “ 0 and so p ■ 1 is an adequate model in that 
case. Still this can be considered a subset of the p • 2 model. There 
does not seem to be a systematic way of determining which, if any, traverses 
have p - 1 short of actually estimating and It is a safe procedure 

to assume p ■ 2 since overfitting will show up the needed simplification. 

In general, it would appear to be a valid and obviously parsimonious step 
to assume at this point that p ■ 2 is a universally applicable part of 
the model. 

At this stage it is apparently appropriate to assume that all three 
parts (p, d and q) of the ARIMA model are as theorized. We, therefore, 
next actually "fit" an ARIMA (2,1,0) model to each traverse. That is, the 
values of and $2 * re estimated from the data using least squares pro- 
cedures. Such procedures are iterative and thus require starting values. 
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Preliminary aatimataa vara obtained from and using the formula given 
previously (page 17). Theaa values ere listed in Table XV. The more 
reliable final estimates are given in Table V. Comparison of these two 
tables shows that in most esses the preliminary and final estisiates ere 
quite close. However, there are enough Instances of large changes that 
reliance should be placed on only the final estimates (Figure 11). An 
additional advantage to obtaining the final estimates is that it allows 
confidence limits to be estimated. This was done and the upper end lower 
95% confidence limits are also given in Table V. These confidence limits 
provide an additional, and very sensitive, test of th appropriateness of 
the postulated ARIMA (2,1,0) model. Having such confidence limits we can 
determine if it is possible to reject the postulated model (as a null 
hypothesis) which would have 0 < < 1 and -1 < $2 < 0. As can be seen 

from Table V all 96 estimated values of fall within the postulated range. 
This should be considered in light of the fact that the po»sible range of 
invertible values is -2 < < 2. More interestingly, only four traverses 

yield values of which are lower than 0.50. This provides a hint that 
the process rates are even more tightly constrained than suggested in the 
original model. Further support for this is given by considering the 
overall mean which is 0.741 with a standard deviation of 0.119. If all 
areas studied are being modified by this process mix at the same rate, 
subject only to normally distributed random fluctuations (an hypothesis to 
be tested shortly), then there is a 95% probability that the true rate is 
in the range 0.717 < $ < 0.765. That is, within the areas studied it 

appears thst these processes (slope wash and creep combined) are changing 
slope form at about three-quarters of the maximum possible rate. Further- 
more, it is highly unlikely (less than 5% chance) to find rates lower than 
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0.50 or greater than 0.98. The highest value actually obaarvad w» 0.98 
(Lavaca, a). The near coincidence of these two upper Halts laply that 
it is highly unlikely to observe non-statlonary slope series; a first 
difference of elevations is adequate and correct, a second difference is 
unnecessary and supercilious. An addltlonsl support for this statement 
is given by the individual upper 951 confidence Units for the estimates 
of jj. In only nine cases does the upper Unit Include 1.0 (the largest 
value is 1.07). Even a value of 1.0 is stationary unless were greater 
than or equal to zero. In no case does the lower 95% confidence limit 
Include zero. Thus, It Is clear that a ^ parameter Is required in the 
nodel. 

Although the hypothesized upper limit for 0^ is 0.0 a fair number of 
traverses ( 31 ) yield values greater than zero. However, In 22 of these 
the value is very close to zero and the lower 95% confidence limit does 
include zero. Thus, taking 0 > # 2 •• th « null hypothesis only nine 
(of 96) traverses allow rejection at the 95% confidence level (about 5 
rejections are to be expected at this confidence level). The nean value 
of all 96 traverses is -0.05 with a standard deviation of 0.126. At the 
95% level of confidence the true population nean nust fall in the range 
- 0.076 < $2 < -0.025 if these are all samples from a single population 
subject only to random fluctuations. Thus, there is no reason to suspect 
that the true expected value varies from the hypothesized range of possible 
values. In each of the eight provinces but one (Blue Ridge) the nean value 
(of 12 traverses) is negative. For the Blue Ridge the value is so small 
(0.011) that it is not significantly different from zero. This would 
seem to imply that within this province (or at least in the three quad- 
rangles studied) creep was not a significant factor in shaping the land form. 
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The presence of some positive values of 4 > ^ i* not necessarily incompat- 
ible with the hypothesized model. If the true value of $2 for a traverse 
is zero and the estimated values are subject to some random (normally 
distributed) error not related to the variability in process rate itself* 
then a positive value is possible. A natural error source in these 
estimates is that resulting from the data collection process itself. A 
considerable number of errors were detected and corrected during the study; 
it is quite conceivable that additional errors remain; although they csn 
be expected tc be small. In some cases errors were not detected until 
estimates of $2 had already been made. When such errors were corrected 

a 


and 4>2 recomputed it was always found that the new value was lower* either 
negative or at least considerably closer to zero. Table XI lists these 
corrections. Such changes are a strong suggestion that sampling errors* 
even minor errors on a single point out of 301, tend to introduce a 
positive bias in the estimate of Because 4^ and 4^ are negatively 

correlated such bias could tend to produce a bias toward smaller values 

A. 

for Such eff^its were also noticed in the corrected series mentioned. 

These effects may also be present in certain traverses. An example of 


such a phenomenon may be given by the Fidelity-c traverse * 0.38, 

/V 

$2 * 0.13); although no error has been detected, the values are anomalous 
and suspicious. 

The correctness of the model is demonstrated quite conclusively by 
the upper 951 confidence limit of 4^. Of the 96 traverses* 27 have an 
upper limit which is negative, indicating that any alternative model is 
highly unlikely fcr at least this proportion of the traverses. Thus, 
the model is sufficient to explain at least 901 of the observations, and 
is uniquely capable of explaining 281 of the observations. 
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Table XI. Effects of errors in traverse elevations upon the final 
estimates of 4>, and Parameters estimated after 
correction of data are also listed for comparison. 


Quadrangle Name 

Traverse 

With 

♦ 

1 

errors 

2 

After correction 

<P d > 

1 2 

Ayer 

c 

.11 

.51 

.77 

.089 

Brandon 

D 

.39 

-.080 

.67 

.059 

Fayetteville 

A 

.50 

.25 

.51 

.25 

Fayetteville 

D 

.67 

.062 

.70 

.033 

Hillsboro 

D 

.57 

-.14 

.57 

-.14 

Lavaca 

B 

.52 

.29 

.53 

.28 

Mena 

A 

.13 

.22 

.76 

-.10 

Mount Mitchell 

A 

.50 

.26 

.75 

.075 

Mount Mitchell 

D 

.52 

.28 

.81 

.017 

Rover 

A 

.38 

.14 

.73 

-.11 

Sherando 

C 

.39 

.39 

.81 

.048 

Strasburg 

A 

.45 

.29 

.75 

.044 
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The preceding analysis shows that when an ARIMA (2,1,0) model is 
assumed the resulting estimates of and $2 within the range 

expected on the bas^ >.■? . omorphic arguments. It remains to be seen 
whether this model with those specific parameter values will actually 
produce a fit to the observed series that is statistically acceptable. 

The goodness of fit of the estimated model can be evaluated by means of 
seven tests; all are based upon the assumption that if the fit is correct 
the residual series will not differ from white noise. This residual 
series is obtained by subtracting values predicted using the model from 
the actual observed series as explained earlier (Chapter II, p. 9 ). 

The results of these tests are described next. 

If the residual series is as expected the mean of the series should 

be aero. A t-test for significant departure from zero is the first test 
(Table XII). As can be seen the null hypothesis (true mean is zero) is 
rejected in only two of the 96 cases (whereas about 5 could be expected 
at the 95% confidence level used). Thus, this test provides no evidence 
to allow rejection of the fit of the model. 

The second test, the portmanteau test, is a check for the tendency 
for significant autocorrelations 1 > occur in the ACF of the residuals. 

If the residuals are white noise the values should not significantly 
differ from zero at any lags. The portmanteau test is based upon the 

Chi-squared statistic and also is reported in Table XII. In general, 

the traverses pass this test, however there are 27 traverses which display 
values exceeding the 95% cu*--off. That this is not a statistical fluke 
is suggested by the fact ti.-i. 10 of the traverses display values exceeding 
the 95.5% cut-off. Thus, there is some indication on the basis of this 
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ORIGINAL PAGE S 
OF POOR QUALITY 

Table XII. Results of first five tests of residuals from the ARIMA 
(2,1,0) model. Tests are explained in the text. An 
asterisk indicates the value exceeds the 95% confidence 
limit . 


i evince 

Quadrangle 

State 

x/st. er. 

Y 2 

X 46 

v^i 

y2p l 

! o 

Interior 

Mammoth 

KY 

a 

.07 

57.3 

-.56 

-.71 


Lev 

Cave 


b 

.26 

43.8 

-.49 

-.67 

.22 




c 

1.10 

35.8 

-.50 

-.65 

.11 




d 

.01 

52.9 

-.52 

-.69 

.23 


Hillsboro 

KY 

a 

.02 

58.1 

-.57 

-.73 

. 2S 




b 

.40 

72.8* 

-.53 

-.71 

.27 




c 

.16 

49.3 

-.52 

-.6?. 

.1*. 




d 

.35 

72.3* 

-.57 

-.75 

.3b- 


Hover 

TA 

a 

.13 

40.4 

-.55 

- .73 

.35- 




b 

.76 

40.6 

-.52 

-.67 

.15 




c 

.44 

39.1 

-.50 

-.67 

.18 




d 

.01 

55.1 

-.50 

-.65 

. 1 3 


Ayer 

MA 

a 

• A A 

68.5* 

-. 70 

- .65 

.0 7 



b 

.25 

67.3* 

-.51 

-.69 

.24 



c 

.27 

67.5* 

-.52 

-.67 

.15 



d 

.05 

44.8 

-.49 

-.65 

. 10 

Kingston 

R1 

a 

1.57 

38.8 

-.53 

-.71 

.32 



b 

.08 

67.5* 

-.49 

-.67 

.24 



c 

.25 

45.6 


-.73 

.3-* 



d 

.60 

84.4* 

-.51 

-.68 

.20 

brand on 

VI 

a 

1.63 

64.5* 

-.53 

- .69 

. i & 



b 

1.78 

55.8 

-.47 

-.66 

.20 



c 

.10 

28.7 

-.50 

-.68 

.21 



d 

00 
«— 1 
• 

— « 

29.5 

-.46 

-.61 

.03 


Uarm 

CA 

a 

1.14 

44.2 

-.54 

-.69 

.20 

Springs 


b 

.14 

42.0 

-.50 

-.65 

.13 



c 

■ .09 

44.1 

-.54 

-.70 

.24 



d 

.78 

46.0 

-.53 

-.65 

.19 

Paterson 

NJ 

a 

.55 

66.8* 

-.53 

-.70 

.73* 



b 

.86 

68.9* 

-.46 

-.62 

.02 



c 

1.56 

33.8 

-.49 

-.66 

.15 



d 

1.03 

41.5 

-.48 

-.64 

.09 
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/rovit.w'S 

Quadrangle 

State 

x/st. er. 

x 2 

X 46 


? 2 Pl 

v 2 p 2 

_edmont 

Washington 

DC 

a 

.02 

60.1 

-.46 

-.63 

.09 


West 


b 

.19 

44.1 

-.47 

-.64 

-.13< 




c 

.83 

50.3 

-.51 

-.67 

.19 




d 

.87 

79.1* 

-.47 

-.64 

.12 

lus Ridge 

Mount 

NC 

a 

.10 

23.2 

-.52 

-.69 

.23 


Mitchell 


b 

.98 

39.9 

-.49 

-.66 

.17 




c 

1.08 

45.9 

-.49 

-.65 

.11 




d 

.33 

50.3 

-.52 

-.69 

.20 


Strasburg 

YA 

a 

1.07 

39.3 

-.51 

-.67 

.16 




L 

.90 

63.3* 

-.51 

-.67 

.14 




c 

.41 

<*S. 2 

-.50 

- .id 

, .i 




d 

.33 

31 .5 

-.48 

-.o4 

. 13 


Pherando 

VA 

a 

.73 

43.3 

-.51 

-.67 

.16 




b 

.72 

54.2 

-.49 

-.63 

.06 




c 

1.30 

45.8 

-.49 

-.65 

.16 




d 

.40 

76.3* 

-.49 

-.66 

.13 

ark 

.renter* 

MO 

a 

1.14 

58.6 

-.51 

-.66 

.15 

riatea^.o 



D 

.36 

46.4 

-.54 

- .7C 

.24 




C 

.13 

47.5 

-.50 

-.67 

.IS 




d 

.26 

58.6 

-.45 

-.64 

.16 


^aint Pa ui 

A K 

a 

.03 

35.6 

-.48 

-.64 

.08 




b 

.41 

47.7 

-.53 

-.70 

.26 




c 

.79 

85.3* 

-.54 

-.71 

.25 




d 

.08 

60.7 

-.54 

-.72 

.29 


> icitiity 

mo 

a 

. 2 C 

56.6 

-.60 

? " 

. j- 




b 

x .80 

36.0 

-.52 

-.72 

.33 




c 

1.32 

100.7* 

-.49 

-.64 

.08 




d 

.16 

31.7 

-.51 

-.67 

.18 


Ouachita Horseshoe Ak a 

Mens . Mm. b 

c 
d 

AK a 
b 
c 
d 


.68 

69.3* 

-.52 

-.69 

.26 

.56 

82.4* 

-.49 

-.64 

.08 

1.16 

68.4* 

-.59 

-.72 

.24 

.47 

61.2 

-.57 

-.73 

.29 

1.46 

41.3 

-.49 

-.65 

.10 

1.53 

44.5 

-.49 

-.66 

.15 

.53 

30.3 

-.50 

-.67 

.13 

.65 

25.4 

-.47 

-.64 

.13 


Mena 
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Province 

Quadrangle 

State 

x/st. er. 

2 

X 46 

V l pj 



luachita 

Lavaca 

AK 

a 

.49 

86.4* 

-.62 

-.78 

.43* 

lit ns. 



b 

1.02 

39.6 

-.51 

-.67 

.18 




c 

1.02 

37.3 

-.45 

-.64 

.14 




d 

1.06 

81.0* 

-.50 

••66 

.17 

'alley and 

Norris 

TN 

a 

.08 

34.3 

-.51 

-.65 

.14 

'■.idge 



b 

.79 

51.4 

-.56 

-.73 

.32 




c 

.25 

58.0 

-.46 

-.60 

-.01* 




d 

.08 

54.6 

-.49 

-.65 

.10 


Alexandria 

PA 

a 

.87 

53.1 

-.50 

-.67 

.18 




b 

.61 

23.5 

-.51 

-.66 

.18 




c 

.60 

41.8 

-.47 

-.65 





d 

.10 

70.3* 

-.42 

-.57 

-.06 s: 


Sau^erties 

NY 

a 

.14 

42.9 

-.44 

-.60 

.01 




b 

.07 

41.3 

-.47 

-.63 

.08 




c 

1.50 

83.3* 

-.47 

-.63 

.04 




d 

.47 

45.6 

-.48 

-.68 

.09 

n-iach . _:r. 

i thaca I'.irst 

NY 

a 

2.77 

61.4 

-.41 

-.56 

0.0* 

. -iceau 



b 

2.78 

35.2 

-.51 

-.69 

.24 




c 

1.54 

62.6 

-.55 

-.71 

.23 




d 

3.53 

38.2 

-.50 

-.69 

.26 


Yavettev -lie 

W\ 

a 

.31 

73.6* 

-.47 

-.64 

.15 




b 

.00 

75.1* 

-.47 

-.65 

.17 




c 

.21 

63.0* 

-.53 

-.72 

.32 

• 



d 

.25 

82.4* 

-.46 

-.64 

.12 


»• hi ewe 11 


a 

. .<6 

7 1 . 7 * 

-.45 

’ • U-» 

• 1 1/ 




b 

.45 

V-T 

OO 

-.48 

-.64 

.08 




c 

. .50 

57.6 

-.47 

-.66 

.19 




d 

.40 

40.7 

-.50 

-.66 

.14 
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test that some structure remains in the pattern of the slope series 
beyond that explained by the slope wash/creep process model. The ACF's 
indicate that this structure must be quite minor. Furthermore, there 
does not appear to be any discernable systematic pattern common to these 
traverses with significant values. For example, one traverse (Brandon-a) 
shows significant values at lags 17 and 29. One would be hard pressed to 
explain such a pattern in any simple manner. Thus, these statistics 
should be taken as a cautionary warning, but in themselves do not suggest 
any systematic alternative model. Of the seven tests of residuals, only 
this one shows any need for caution. 

The third test is based upon the fact that when the ACF of the first 
difference of a white noise series is calculated the values at all lags 
will be zero except for lag one, where the value -0.50 will be observed. 
Thus, the ACF of each traverse can be tested to see if that lag value 
differs significantly from -0.50. An approximate confidence interval 
can be constructed using as standard error the value obtained under the 
null hypothesis that the first difference series Itself is white noise. 
That standard error is 0.06, thus an approximate 95% confidence interval 
is -0.38 to -0.62. As can be seen in Table XII, none of the 96 traverses 
yield values outside this range. Thus, this test gives no reason to 
suspect the goodness of fit of the estimated models. 

The fourth test arises because the second difference of a white 
noise series will also display a specific value at lag one, namely -0.67. 
An analogous test is constructed and in this case the standard error is 
identical, 0.06. Thus, a 95% confidence interval is -0.55 to -0.79. 


Again, it was found that no traverses yield values outside of this range 
(Table XII). The overall goodness of fit is further reinforced. 

The fifth test is also based upon the expected value of a second 
difference of white noise. In this case it can be shown that, at lag 
two, a value of 0.17 should appear. All higher lags have an expectation 
of zero. The approximate standard error in this case is 0.08 giving a 
95% confidence Interval of 0.01 to 0.33. This test resulted in signifi- 
cant values in 10 of the 96 traverses (Table XII). There does not appear 
to be any tendency for the values to be too high or too low as six values 
exceeded the upper limit and four fell below the lower limit. This 
number of significant values does not appear to be beyond the acceptable 
limit for a 95% confidence level which is based upon an approximation of 
the true standard error. Thus, with 86 of the 96 traverses easily passing 
this approximate test there appears to be no reason to suspect the good- 
ness of fit of the model on the basis of this test. 

Summarizing the five tests we have three which reject fewer than 
would be expected at the 95% level (l.e. zero as opposed to 5), one which 
rejects about what would be expected, and one which rejects more than 
would be expected. Combining the results of applying all five tests, 
giving a total of 480 tests, we see there were 37 rejections for an overall 
average of 8%. This is so close to the nominal 95% level that there appears 
to be no good reason to doubt the overall ability of this model to provide 
adequate goodness of fit to arbitrary traverses. Inspection of Table XII 
will show that there is no systematic tendency for a particular province 
or quadrangle or traverse orientation to fail to fit the model. Thus, 
it seems reasonable to conclude that the results obtained are precisely 
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the kind of random variability one would expect at the 951 confidence 
level If the null hypothesis were indeed true. The goodness of fit of 
the models is as desired. 

Tests six and seven are both based upon the partial autocorrelation 
function. Again they are based upon the assumption that the residual 
series is not significantly different from white noise. Test six makes 
use of the fact that the PACF of the first difference of a white noise 
series will display a value of -1.0 at lag zero and equal -l/(k+l) at 
successive lags, k. An approximate standard error (based on the assumption 
that that series itself is white noise) is given by 1/ N where N is the 
length of the series. Two standard error (approximately 95Z confidence 
level) limits were computed for lags 1 through 7. Beyond the 7th lag the 
values are not significantly different from zero. Lags for which the partial 
autocorrelation fell outside of these limits are listed in Table XIII. 

Since each of the seven lags yields a test, a total of 672 tests are 
provided. As can be seen only 42 lags or about 6% are significant. This 
is quite reasonable for a test designed to reject approximately 5Z of the 
results when the null hypothesis is true. Test six does not provide any 
evidence sufficient to reject the goodness of fit of the models in general. 
More importantly, no single test yields more than two significant values. 
There does not appear to be any reason to reject the goodness of fit of 
the model in any specific instance. 

Test seven is available because the PACF of the second difference 
of a white noise series will decay according to the formula -2/(k+2) at 
each lag k. Again approximate 95% confidence limits were set up and values 
of lags one to seven outside these limits are listed in Table XIII. Only 
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Table XIII. Lags (out of first seven) which are significant InPACFof first 
difference (Vl) and second difference (V?) of white noise series. 
Numbers should be read as: 3 means the third leg is significant, not 

that three lags are significant, zero means no lags are significant. 


Traverse 


A 


B 


C D 



Province 

! 

Quadrangle 

V 1 

7 2 

V 1 

V 2 

V 1 

V 2 

7 1 


1. 

Interior 

Mammoth Cave 

3 

0 

0 

6 

0 

0 

6 

6 


Lowlands 

Hillsboro 

0 

0 

0 

7 

6 

6 

2 

0 



Rover 

2 

5 

0 

0 

0 

0 

0 

0 

2. 

New England 

Ayer 

3 

3.7 

0 

0 

7 

6 

0 

0 



Kingston 

0 

0 

5 

5,7 

2 

7 

0 

0 

j 


Brandon 

3 

0 

0 

0 

0 

0 

0 

0 

3. 

Piedmont 

Warm Springs 

0 

0 

5 

0 

0 

0 

0 

0 



Paterson 

6 

5 

0 

6 

0 

0 

0 

0 



Washington West 

0 

0 

0 

0 

5 

0 

7 

0 

4. 

Blue Ridge 

Mt. Mitchell 

0 

0 

0 

0 

0 

0 

0 

0 

I 


Strasburg 

0 

0 

0 

6 

0 

0 

0 

0 

1 


Sherando 

0 

0 

0 

6 

0 

0 

3 

3 

5. 

Ozark Plateaus 

I ronton 

5 

7 

7 

0 

0 

0 

4 

0 



St. Paul 

0 

7 

4 

0 

0 

0 

0 

0 



Fidelity 

0 

1 

2,6 

6 

0 

0 

0 

0 

6. 

Ouachita Mtns. 

Horseshoe Mtn. 

0 

0 

0 

0 

3.5 

5 

3,5 

0 



Mena 

0 

0 

0 

0 

3 

0 

0 

0 



Lavaca 

2 

3 

0 

0 

0 

5,6 

5 

7 

7. 

Valley and 

Norris 

0 

5,6,7 

0 

0 

5 

5,6 

0 

0 


Ridge 

Alexandria 

0 

0 

0 

7 

0 

4,6 

2.4 

0 



Saugerties 

0 

0 

o 

0 

5 

0 

4 

0 

£. 

Appalachian 

Ithaca West 

0 

0 

0 

0 

3 

o 

0 

0 


Plateaus 

Fayetteville 

6 

0 

0 

0 

2 

0 

0 

0 



Whitwell 

0 

! 

0 

0 

o 

4 

0 

0 

6 
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34 of 672 value* fall outsidt the expected range. Since about 5X era 
significant and about 5% would be expected at this confidence limit there 
is no evidence to reject the goodness of fit of the models in general. 
Again, there is also no specific evidence to reject the goodness of fit 
in any particular case. It appears that all traverses are modelled 
adequately using the ARIMA (2,1,0) model and th* estimated parameters. 

Because tests six and seven sctually comprise 1344 separate tests 
a further examination of the results can be made which sheds light on the 
actual results. To motivate this test we note that one traverse (Norrls-a) 
actually yielded three significant lags; is this enough to suspect the 
goodness of fit of the Norris-d model? After all this is a 43% rejection 
rate. First note that any rejection, even one, would exceed the 5% 
rejection rate expected (giving a 14% rate). And yet in 672 tests we 
can expect many to reject one lag. Thus, it is the overall pattern 
which is most important. Part of that pattern is given by the number of 
times no lags, one lag, two lags, three lsgs, etc. were actually rejected. 
Under the null hypothesis and at the confidence level chosen *ach test 
may reject or not with a probability of 0.05. An idea of how many times 
0, 1,2,3, etc. rejections can be expected even if the null hypothesis is 
correct can be obtained by application of the binomial theorem. Ue take 
each traverse ACF to be a binomial experiment consisting of seven trails. 

A total of 192 such experiments have been performed. With N“192, n*7, 
p“0.05 and q»0.95 we will expect the following results: 
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number of lags number number 

slgnif leant expected observed 

0 134 132 

1 49 50 

2 8 9 

3 1 1 

192 192 


As can be seen there Is no Indication of any pattern in the observed 
results. They can be considered independent realizations of a random 
variable (exceeds 0.95 or limit or not) which has p*0.05. This test 
provides further firm evidence that the model fits the observed data 
well. 

Having demonstrated that the AR1MA (2,1,0) model is reasonable and 
that the estimated parameters are in general compatible with the expected 
range, it is now appropriate to summarize the estimates of these parameter 
values. In Table XIV are the means of the four traverses in each quadrangle. 
Also listed are the mean values for each province and the grand mean 
computed from all 96 traverses. As can be seen, there is a considerable 
range in the mean values for provinces (0.64 < < 0.84, -0.12 < ^ < 0.01) 

which in general is considerably larger than the variance within provinces. 

A 

This is especially noticeable in the case of the 

There is obviously the suggestion that there are significant differences 
from province to province in the values of and possibly of $ 2 * One must 


Table XIV. Summary of parameter values for each quadrangle 


PHI 1 


QUAD 


Province 

1 

2 

3 

MEAN 


Interior Lowlands 

0.75 

0.77 

0.74 

0.75 


New Englend 

0.78 

0.76 

0.82 

0.79 


Piedmont 

0.86 

0.80 

0.80 

0.82 


Blue Ridge 

0.84 

0.64 

0.68 

0.72 


Orark Plateaus 

0.82 

0.71 

0.72 

0.75 


Ouachita Mtns. 

0.69 

0.61 

0.62 

0.64 


Valley and Ridge 

0.75 

0.76 

0.66 

0.72 


Appalachian Plateaus 

0.82 

0.71 

0.72 

0.75 


QUAD MEANS 

0.79 

0.72 

0.72 



GRAND MEAN 




0.74 



-.04 
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•Iso consider the possibility thst within provinces there nsy be significant 
differences from qusdrsngle to qusdrsngle. Xt is also of soas Interest to 
svsluste whether there is • significant tendency for one traverse orienta- 
tion to vary from another orientation in terns of these parsswter values. 
These questions have been examined by means of an analysis of variance 
performed on each parameter separately. 

The analysis structure is as follows. There are three factors whose 
effects upon the parameters is of interest: province, quadrangle (l.a. 

ares within province), and traverse (i.e. orientation). Province is 
measured at eight levels, quadrangle at three in each province, and 
traverse at four within each quadrangle. The analysis is a random 
effects model for all three factors; if it were redone, new levels would 
almost surely be included. The quadrangle factor muat be considered 
nested within the province factor since it is not possible to find the 
same quadrangle in some other province. Although traverse could be 
considered to be nested within quadrangle for this study the orientation 
of the traverse was the phenomena of interest. In that case a traverse 
of each orientation can be observed in each province at each quadrangle. 
Thus, the factor traverse was considered to be crossed with province, so 
that in addition to province, quadrangle nested in province, and traverse 
another source of variation whose effect was computed was the province 
x traverse interaction. The error term was estimated by the quadrangle 
by traverse interaction nested within province which had 48 degrees of 
freedom in each test. This was used to test both the province by traverse 
interaction and the traverse main effect. A second error term, the 
quadrangle nested in province effect, was used to test the main effect 
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of province unless it was found to be not significant, in which case 
the other error term was also used for this test since with more degrees 
of freedom it could be expected to provide a better estimate of the 
true error term. The computed analysis of variance tables are given 
(Tables XV and XVI). 

In the case of only one source of variation was found to be 

significant. That is, the main effect of province. All other sources 
of variation yielded remarkably consistent mean squares which were nearly 
identical to that estimated from the "error" term. 

The analysis of variance for $2 shows that the quadrangle effect is 
significant. The mean square due to province is actually slightly large' 
but because its expected mean square includes a contribution from the 
quadrangle effect that term was used as the error term in computation of 
the F-ratio. It was thus found that the province factor is not significant. 
Differences between provinces are no greater than can be found within a 
single province in going from one area to another. Note that, according 
to the hypothesized model, the coefficient phi 2 is equal to the ne<^t.ive 
of the rate of creep, thus these results apply to the creep rate coefficients 
equally as well. 

Since the hypothesized model is so well supported by the available 
data, it is reasonable to explicitly break out the slope wash rate 
coefficient and study its variability from area to area. Recall that phi 
1 is the total rate equal to the sum of the creep and slope wash rates. 

Thus, the sum phi 1 plus phi 2 is an estimate of the slope wash rate. 

This rate was calculated for the 96 traverses And an analysis of variance 
run using the same design as previously. Table XVII shows the results. 
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Table XV . Analysis of variance for PHI 1. An asterisk ("*") indicates signif- 
icance at the 95Z confidence level , "N.S." means not significant. 


Source of 
Variation 


Degrees of 
Freedom 

Sums of 
Squares 

Mean 

Square 

Error 

Term 

F 

Ratio 

Signif- 

icance 

Province 

A 

7 

.241 

.034 

CB(A) 

2.43 

* 

Quad (A) 

C(A) 

16 

.201 

.013 

CB(A) 

.93 

N.S. 

Traverse 

B 

3 

.027 

.009 

CB(A) 

.64 

N.S. 


AB 

21 

.227 

.011 

CB(A) 

.79 

N.S. 

"Error” 

CB(A) 

48 

.649 

.014 




Total 


95 

1.345 

.014 





Table XVI. 

Analysis of variance for PHI 2, 
rate coefficient. Abbreviations 

equal to the negative of the 
as in previous table. 

creep 

Source of 
Variation 


Degrees of 
Freedom 

Sums of 
Squares 

Mean 

Square 

Error 

Term 

F 

Ratio 

Signif- 

icance 

Province 

A 

7 

.244 

.035 

C(A) 

1.17 

N.S. 

Quad (A) 

C(A) 

16 

.485 

.030 

CB(A) 

2.73 

* 

Traverse 

B 

3 

.010 

.003 

CB(A) 

.27 

N.S. 


AB 

21 

.247 

.012 

CB(A) 

1.09 

N.S. 

"Error" 

CB(A) 

48 

.517 

.011 




Total 


95 

1.504 

.016 
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As with creep, the quadrangle tens Is found to be significant indicating 
that creep varies from area to area within a single physiographic pro- 
vince. As with creep, any apparant differences between provinces can be 
attributed to the particular sites selected for comparison. Again, 
travese orientation has no effect on its own, and the interaction term 
(traverse by province) also shows that no differences in slope wash rate 
estimates results even when comparing one traverse direction in one pro- 
vince to a different traverse direction in some other province. 

To illustrate the differences in the values of phi 1 (overall process 
rates) that can be seen between provinces and to determine which, if any, 
provinces could be considered to behave similarly, a multiple comparison 
test utilizing the Least Significant Difference method (Figure 15 and 
Table XVIII). As can be seen in the table, at most three distinct group- 
ings of overall process rates could be defined. As can be seen from the 
overlap in group assignments, no distinct separation is possible and it 
might appear just as reasonable to consider the variation to be a 
continuoum. A more conservative comparison procedure, Tukey’s "Honestly 
Significant Difference" technique suggests that at most two groupings 
could be defined (Table XIX). It is important to note that neither method 
would allow all provinces to be assigned to a single group. However, 
there is one (extremely) conservative multiple comparison method, namely 
Sheffe's Test, which considers all possible comparisons, that does group 
all provinces together. This is not consistent with the ANOVA results 
and one of the other groupings, probably the middle one, Tukey’s HSD, 
is to be preferred. Following that method, we consider the Ouachita 
Mountains to have a low process rate and the remaining provinces to have a 
(comparatively) high process rate. 


Table XVIII. Least significant difference test on 
parameter 4>j« 


Level 

Multiple Comparisons 
Mean 

on Factor PROVINCE 
Semple Sire 

Separation 

6 

.64 

12 

a 

4 

.72 

12 

ab 

7 

.72 

12 

ab 

8 

.75 

12 

be 

5 

.75 

12 

be 

1 

.75 

12 

be 

2 

.79 

12 

be 

3 

.82 

12 

c 


Least Significant Difference 

Error mean square ■ .014 
Degrees of freedom * 48 
Alpha level ■ .05 

Table value from Student's t • 2.01 
LSD value - 9.709 
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Table XIX. Tukey's Honestly Significant Difference test 
on parameter 4) . 


Level 

Multiple Comparisons 
Mean 

on Factor PROVINCE 
Sample Size 

Separation 

6 

.64 

12 

a 

4 

CM 

• 

12 

ab 

7 

.72 

12 

ab 

8 

.75 

12 

ab 

5 

.75 

12 

ab 

1 

.75 

12 

ab 

2 

.79 

12 

ab 

3 

.82 

12 

b 


Tukey's HSD 

Error mean square » 0.014 
Degrees of freedom * 48 
Alpha level ■ .05 

Table value from Studentized range ■ 4.48 
HSD value - .153 
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v. INTERPRETATIONS 


Because of the speculative nature of this study there are a number 
of ite is that should ba described separate from the results per se . These 
are separated from the co nclusions section because they are isore pre- 
cisely considered to be geemorpnlc interpretations . Conclusions pertinent to 
the geobotanJ. apnlicatf .»ns intended in this study — which are, in large 
part, based upon these gc'morphic Interpretations — will be discussed in 
a succeeding chanter. 

When '..as .vati . at the scale of this study it is definite that there 
are a-, jure^ble relations between adjacent slopes in a traverse of the 
.a.tdfo”;n. Such relations are appropriately described by the postulated 
ARINA (2,1,0) model. Using this model good estimates of the rate of 
landform modification by the two processes - slope wash and creep, sensu 
lato - may be obtained. In particular, -$ 2 , is an estimate of the rate 
of creep within the area. Slope wash rate is given by . Thus $1 

is an estimate of the overall rate of landform modification. Because of 
the considerable variability in these rates which can be observed to occur 
even within a small area (the size of a quadrangle), it is best to use the 
averages of several traverses when comparing regions. In general, it 
appears that the rate of landform modification by slope wash greatly 
exceeds the rate of creep. Even in a single traverse the rate of slope 
wash is usually twice that of creep. That the rate of slope wash 
exceeds that of creep agrees with the findings of Carson and Kirkby 
(1972, p .tlf). No traverse in the entire group studied yielded estimates 
of creep as large as that of slope wash. Variability in the rates of 
the two processes is apparently about the same. 
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The analysis of variance has shown that real differences do exist In 
the degree of relation between adjacent slopes In the various regions 
studied. The overall rate of landfom modification (as measured by phi 1) 
varies from province to province. There Is an Indication that thare is a 
continuous variation in such rates. However, the Tukey HSD test Implies 
that at least two distinct levels of rates are found. A reasonable 
interpretation is that the Ouachita Mountains province Is modified at a 
distinctly lower rate than the remaining provinces. Such differences 
apply to the overall rates but cannot be ascribed to either slope wash 
or creep specifically, only to their combined activity. 

Analysis of phi 2 (creep) and the sum phi 1 plus phi 2 (slope wash) 
shows that the variability of these rates is at a finer scale than the 
province level. Significant differences in these rates cannot be demon- 
strated between provinces but do exist from quadrangle to quadrangle within 
provinces. Such differences are much larger than the province to province 
variations. Thus, there are regional differences in the rates of both 
creep and slope wash which do not show any systematic relation to province. 
This could very well reflect local variations in controlling factors such as 
lithology or structure as well as the variations in climate. Interesting- 
ly, the three analyses together imply that the overall process rates are 
quite constant within any province and since both of the component rates 
do vary significantly within the province, it seems that decreases in one 
rate within a province are made up by increases in the other. As will be 
seen in the final chapter there appears to be a negative relation between 
slope wash and creep rates, thus it would appear that overall process rate 
may be a significant basic (and perhaps underlying) feature of each physio- 
graphic province. 
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Importantly, it la clear that the degree of relation between adjacent 
slopes estimated (i.e. the rates) do not depend at all upon the orienta- 
tion of the traverse being considered. This is true for phi 1, phi 2, 
slope wash and creep. Such a finding means that detailed study of rates 
within an area the size of a quadrangle is not required. Even one traverse 
is sufficient. Perhaps study of a single traverse in each of several ad- 
jacent quadrangles would yield better estimates of the regional rates. 

Finally an important observation arising from these studies is a 
distinction between form and process. In particular, it is important to 
recognize that the techniques discussed yield estimates of the rates of 
processes responsible for forming the slopes. These are always past 
rates which may or may not correspond to present rates. Rates may 
change considerably faster than the slope forms can respond. An example 
may be given by the Brandon, Vermont quadrangle in the New England physio- 
graphic province. This is an area of sharp relief significantly carved 
by continental glaciation during Late Wisconsinan times. These forms are 
only partially modifi(! by subsequent processes. On the other hand, phi 1, 
the estimate of overall process rate, is 0.82 for that quadrangle. This 
is a very high value, exceeded in only two other quadrangles of the 24 
studied. Such a value could very well reflect higher process rates shaping 
the landform under (perhaps) periglacial conditions. 
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VI. CONCLUSIONS 

Process rates effecting the regolith cen be estimated frost eveileble 
topogrephic steps. In general, mean values will be about 0.75 + 0.15 for 
any area and all values within the eras will be within + .12 of one 
another. Such rates give an estimate of how rapidly the slope forms are 
being modified by slope wssh and creep. This, in turn, should be related 
to the rate of weathering of the underlying regolith and bedrock. If 
this bedrock contains ore bodies the strength of their signal In the 
regolith can be expected to be Influenced strongly by the rate of which 
that regolith is developing. Thus, we can expect that our ability to 
use geobotanical exploration will be related to the measure of overall 
process rate given by phi 1. 

Although overall process rates remain at a consistent level within 
a given physiographic province (l.e. there is no significant quadrangle 
effect on phi 1) the process mix can still vary significantly from area 
to area within that province. Thus, both slope wash and creep vary 
significantly from area to area within a given province. The particular 
mix that occurs in an area can be expected to strongly influence the 
strength of a geobotanical signal of an oze body. In particular, slope 
wash tends to carry material far dovnslope and thus spread it over a large 
area. In contrast materials are moved slowly and only short distances 
by creep. Thus, in areas in which creel is relatively strong the geo- 
botanlcel signal of an ore body should also be relatively strong. Where 
the slope wash effects are especially large the rate of regolith removal 
may be to great to allow significant geobotanical signals of ore bodies 
to build up. Thus, such investigations should be emphasised in regions 
with a large (in a negative sense) value of phi 2. 


no 


Because of phi 2 term is so snail, and in many casts is not signif- 
icantly different from zero, it is apparent that the AK(1) portion of the 
model dominates (in a slope model). Thus, the arguments of Craig (1982) 
probably give a very close approximate measure of the scale of the land- 
form in a given area. It has been shown that this terrain scale also 
sffects the properties of remote sensing data. That the phi 1 term 
dominates in defining landform scale seems quite reasonable in that it 
is a measure of the overall rate at which processes shape that fora. 

In general, the more active the processes the larger an area the ore body 
signal can be expected to be spread over. Thus, the spsclng of camples 
in geobotanical investigations should increase as the value of phi 1 
Increases. Estimates of the spacing using the value of phi 1 can be 
computed from the formula in Craig (1982). 
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VII. FUTURE WORK 

On the basis of the concepts and results discussed above six major 
areas of additional work can be defined. At present so little work has 
been done relating form and process in a quantitative manner that our 
understanding is perhaps just at a threshold. Whether significant pro- 
gress can be expected in the future depends upon how clearly the path of 
future investigations is defined. Some gaps in our knowledge are listed 
below. 

Ensuring the accuracy of the data collection method is a critical 
need. It has been shown that a large number of errors naturally arise 
using the methods now employed. These methods include interpolation 
between contour lines, physically recording data and manual keying it 
into the computer. Each of these steps introduces errors which are difficult 
and time-consuming to detect and correct; and of course this necessarily 
introduces a new step in the procedure. An example of the enormous 
effect of a single error will illustrate the need to avoid such errors 
and minimize their probability of occurrence. Traverse C of the Ayer 
quadrangle contained an error at point 52. Table XX summarizes some 
parameters computed from the data before and after the correction was 
made. It is clear that such an error is sufficient to make the results 
useless. It is fortunate that such errors have an effect which is dramatic 
and hence makes them obvious. 

Three techniques to avoid such errors have been considered and are 
presently undergoing evaluation. Each is based on the idea that the 
data should go directly into the computer in one step rather than a 
multi-step operation. The first method would make use of the digital 
terrain models available from the United States Geological Survey (U.S.G.S.) 


Table XX. An example of the effects of a single error 
upon the estimates of various statistics. 
Traverse C, Ayer HA. 


Parameter Error Corrected 


Rho (1) 

0.303 

0.818 

Rho (2) 

0.324 

0.700 

Estimate of 

0.82 

0.74 

Phi (1) 

Estimate of 

-0.09 

0.10 

Phi (2) 

Final Phi (1) 

0.11 

0.77 

Final Phi (2) 

0.51 

0.09 


Chi-squared 
for ACF 


23.7 


67.5 


for some quadrangles. Such models are contained on magnetic tapes and 
contain enough data, in grid form, to construct the topographic maps. 

This data can be extracted directly and input to the autocorrelation 
programs. Unfortunately, only a small percentage of the total area of 
the U.S. has yet been so mapped. And, of course, virtually nowhere else 
in the world will such maps be available in the foreseeable future. Con- 
sequently, other error-free methods must still be available. 

A second method is to use a digitizer to follow the same traverse. 

This time instead of equal intervals we digitize the points at which 
the contour lines cross and indicate whether we are going uphill or down. 
From this information the computer can be used to reconstruct the cross- 
section and exact linear interpolation can be used. Even more sophisti- 
cated interpolation is easily achieved; presently cubic interpolation 
appears to be most desireable. 

Another procedure being explored is to interpolate mentally as before 
but to input the resulting data directly using speech recognition techni- 
ques. Although slower and probably less accurate than digitizing, it is 
still much faster and more precise than the old manual method. 

Besides the means of data input, data collection can also be improved 
in several ways. It would be useful to determine the effects of changing 
the sample spacing within a traverse. Perhaps the same level of preci- 
sion could be achieved using 2 mm rather than 1 mm spacing. There is 
some indication that changing the sample spacing will have a systematic 
and predictable effect upon the parameter estimates. This is true for a 
pure ARIMA (l,d,0) model (Craig, 1982). Whether it can also be demonstrated 
for an ARIMA (2,d,0) model needs to be examined. 


. 
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Other changes in data collection could include taking fewer traverses 
per quadrangle, since one is apparently enough and orientation does not 
introduce a bias. Whether replicates in adjacent quadrangles are de- 
sireable should be examined. The large effects of errors makes it 
desireable to introduce more tests for errors in a systematic fashion. 
Since errors in effect add a 'noise' series to the true data a more 
complete analysis of the theoretical effects of added noise upon the 
model structure would be of use. 

A second major area of concern in future studies should be the 
extent to which ‘■hese results and in particular the proposed ARIMA (2,1,0) 
process/form model can be extended to other area. At a minimum samples 
should be taken from each physiographic province in the U.S. Prefer- 
ably, estimates of process rates should be available from every physio- 
graphic section in the U.S. Such a study is, in fact, presently underway 
(NASA Grant NAG 5-166 to KSU). These data should make it clear whether 
mass wasting processes (mudflows, rockfalls, landslides, slumps, e*-c.) 
add a significant new structure to the autocorrelational properties of 
slopes. Some samples will be taken from areas known to be dominated by 
such processes. 

In general, it would be wise to apply this method of estimating 
process rates to areas where previous studies have documented the actual 
process mix based upon field observations. A number of such studies, 
spread throughout the world, have been documented in the literature. 

The two methods of estimation should be compared directly. 

Such studies are relatively rare and tend to concentrate on one 
process rather than a total evaluation. Thus , the most critical test of 
the postulated model will probably only come if field examinations are 
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specifically designed to evaluate the accuracy of its predictions in the 
field. Perhaps the simplest of such tests would arise when field work 
occurs in areas of extremes of the processes. Examples which could be 
appropriate Include: badlands areas, areas of dunes, trop.'.cal regions 

and the dry valleys of the Antarctic. Estimation of process rates in 
the field is a difficult and time-consuming task. Because of the seasonal 
and yearly variations in the true rates such studies would ideally extend 
over a number of field seasons. The techniques of such studies are 
fairly well worked out (Goudie, 1981). It will take a significant 
committment of time and resources to achieve this objective. 

We can expect that a number of studies of geomorphic interest can 
be pursued with such process-rate data available. As an example, consider 
the question: are the rates of surface wash and creep related? The 

results of the analysis of variance tests had suggested that an increase 
in one should be accompanied by a decrease in the other so that the sum 
for a province remains constant. This suggests the following, let k be 
the province sum, b the creep rate and a the rate of slope wash, then: 

k * a + b 
a - k + (— 1 ) *b 

We have a classic regression model relating the two rates in which the 
intercept equals the province sum and the slope of the line is minus one. 

A plot of the actual data (Figure 16) shows the inverse relation as 
expected. However, the estimated slope is not -1, rather it is -.56 and 
the intercept does not equal the mean overall rate .74, rather it is .43 
(Table XXI). Perhaps this discrepancy arises because all 96 traverses of 
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SLOPE NOSH AND CREEP RELATION 



SLOPE WASH 



Table XXI. 

Source of 
Variation 

Linear 

Quadratic 

Residual 


Analysis of variance table of polynomial regression 
of creep rate (dependent variable) versus slope wash 
rate (Independent variable). 


Degrees of 
Freedom 

Sums of 
Squares 

Mean 

Square 

F- 

Ratio 

1 

.442 

. 442 

43.08 

1 

.108 

.108 

10.52 

93 

.954 

.010 



TOTAL 


95 


1.504 


* 
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data arc lumped together in one regreaslon. The ANOVA reaulte would 
auggeat that the relation would ahow up when conaiderlng the data of 
individual provinces. It should be examined in more detail. Surprisingly, 
the linear regression is significantly Improved upon by a quadratic one 
(Table XXI) . This seems to suggest that creep reaches a maximum when 
surface wash takes a middle value and decreases at lower or higher 
values. This is conceivable in that high values of surface wash might 
remove material before creep can move it and low values of surface wash 
(presumably in an area of low rainfall or one lacking regollth) would be 
accompanied (in a dry area) by lack of moisture to lubricate and enhance 
the creep process. Whether the relation is truely quadratic cannot be 
demonstrated until more data from areas of low slope wash rates are 
available. If the quadratic result is shown correct on the basis of 
additional data, one can expect at least one ANOVA result to change also. 
It is clear that many interesting geomorphic studies are possible with 
such data. 

Additional geomorphic questions which can be examined with such 
data should emphasize the relation of the process rates to variables 
which could be responsible for them. Such variables include: 


i) temperature 
ii) precipitation 
ill) elevation 
iv) thickness of regollth 
v) erodibillty of bedrock 


vl) vegetation cover 
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Certain of that* data art aaally collect ad. Othara require much sore 
work. Especially the latter three would probably require field work to 
catenate well. It Is natural to expect these process rates to be re- 
lated to the rate of denudation in an area. Technically, the rates 
estimated using the ARIMA (2,1,0) model are the rates these processes 
acted at during the time that the landform was shaped to Its present 
form. In spite of the technical difference a relation between process 
rates and denudation rates would not be surprising to see and would 
certainly be useful to know of. Another question of considerable Interest 
is whether process rate estimates will differ depending on whether the 
slope processes are dominantly transport limited or weathering limited. 

The autocorrelations! structure of these two classes of slopes should 
be compared and contrasted. 

It is now clear that slope traverses follow a specific model of 
autocorrelation. This means that the angle of a slope is related to 
thit of adjacent slopes. Because slope angle has a significant effect 
upon the distribution of reflected light from that slope, it can be 
expected that sequences of reflected light measurements along a traverse 
of the ground (such as represented by LANDSAT data) will also show an 
autocorrelation structure. That such a structure actually exists has 
been demonstrated (Craig, 1979; Craig and Labovitz, 1980). It also 
appears that it is related to properties of the terrain below (Craig, 
1981). It remains to be demonstrated that there is a direct link between 
the two data structures. This would be appropriately studied by means 
of transfer function theory (Box and Jenkins, 1970). Such a study remains 
to be done but is now appropriate. 
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